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PREFACE 
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Supply  Center  (DESC),  Dayton,  Ohio,  serves  as  the  contracting  agency.  The  ASC 
program  manager  was  Doug  Amlin.  The  ERIM  program  manager  was  Michael  T. 
Eismann.  The  authors  of  this  report  are  Michael  T.  Eismann,  Stephen  D.  Ingle,  and  M. 
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1.0  INTRODUCTION 


The  Image  Based  Sensor  Model  (IBSM)  is  a  modular  set  of  numerical  tools  for 
designing,  evaluating,  and  modeling  electro-optical  and  infrared  (EO/IR)  imaging  sensors. 
The  primary  motivation  which  led  to  the  development  of  IBSM  was  the  need  for  a  model 
which  (a)  could  produce  simulated  sensor  imagery  in  addition  to  sensor  performance 
metrics  to  better  characterize  the  imaging  performance  of  a  sensor  system,  and  (b)  provides 
the  flexibility  to  evaluate  and  compare  sensors  and  imaging  configurations  with  differing 
characteristics  without  rewriting  computer  code. 

To  satisfy  the  first  need,  IBSM  is  primarily  geared  to  image  simulation;  that  is,  it 
can  generate  renditions  of  high  fidelity  input  scenes  which  are  representative  of  actual 
sensor  imagery,  including  a  wide  range  of  atmospheric,  sensor,  data  link,  and  processing 
effects.  In  addition  to  this  image-in,  image-out  mode  of  operation,  the  model  also  contains 
tools  for  parametrically  computing  sensor  performance  metrics  such  as  signal-to-noise  ratio 
(SNR),  ground  sample  distance  (GSD),  ground  resolved  distance  (GRD),  and  image 
quality  as  defined  by  the  National  Image  Interpretability  Rating  Scale  (NIIRS).  The 
combination  of  these  two  operating  modes  provides  the  user  the  means  for  evaluating  the 
level  of  synergy  between  subjective  imaging  performance  to  quantitative  sensor 
performance  metrics. 

To  satisfy  the  second  need,  IBSM  was  designed  to  be  completely  modular  in  that 
each  effect  modeled  is  done  so  as  a  self-contained  function.  This  provides  the  user  with  a 
set  of  modeling  building  blocks  which  can  be  tied  together  in  a  wide  variety  of  ways  to 
easily  model  almost  any  sensor  design  or  imaging  configuration.  Each  module  contains  its 
own  graphical  user  interface  (GUI)  for  simple  entry  of  the  defining  module  parameters. 

The  modular  design  philosophy  also  has  the  advantage  that  it  easily  accommodates 
upgrades  to  extend  it  to  sensing  and  processing  effects  not  currently  covered.  This  simply 
involves  the  generation  of  a  new  building  block  as  opposed  to  a  full  model  revision.  As  an 
example,  the  current  model  does  not  contain  tools  for  synthetic  scene  generation,  but  such  a 
capability  could  be  added  without  impacting  the  remainder  of  the  model. 

This  report  overviews  the  IBSM  model,  including  a  brief  overview  (Section  2), 
detailed  description  (Section  3),  usage  overview  (Section  4),  example  models  (Section  5), 
and  software  overview  (Section  6). 
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2.0  MODEL  OVERVIEW 


2.1  MODELING  APPROACH 

IBSM  consists  of  an  extensive  set  of  discrete  modeling  tools  which  can  be  tied 
together  to  evaluate  the  performance  of  a  wide  variety  of  sensors.  IBSM  was  primarily 
intended  to  be  used  in  an  image  simulation  mode  as  depicted  in  Figure  2-1.  Such  a 
simulation  begins  with  a  high  fidelity  input  image  representing  the  target  and/or  background 
scene  to  be  simulated.  The  image  is  then  degraded  through  a  series  of  operations 
characterizing  the  intervening  atmosphere,  the  sensor  optics,  the  electro-optical  detector,  the 
data  link,  and  subsequent  signal  processing.  Various  utilities  for  image  format  conversion, 
image/spatial  spectrum  conversion,  metric  computation,  and  image  display  also  exist. 

A  parametric  modeling  capability  is  also  supported  by  IBSM,  specifically  for 
computing  image  quality.  This  is  depicted  in  Figure  2-2.  As  opposed  to  the  image 
simulation  mode,  parametric  modeling  does  not  explicitly  involve  an  image  chain.  Rather, 
sensor  metrics  such  as  the  modulation  transfer  function  (MTF),  signal-to-noise  ratio 
(SNR),  ground  resolved  distance  (GRD),  and  image  quality  (NIIRS)  are  computed  directly 
from  sensor  parameters.  This  computation  can  be  performed  iteratively  while  altering  a 
specified  parameter  within  the  loop  to  investigate  performance  as  a  function  of  an  arbitrary 
independent  variable. 

In  general,  the  model  incorporates  multispectral  imagery  (MSI).  In  the  image 
simulation  mode,  MSI  are  explicitly  carried  through  the  image  chain.  AdditionaUy,  the 
center  wavelengths  and  spectral  bandwidths  of  each  band  are  carried  through  the  chain  such 
that  wavelength  dependent  modules  operate  differently  on  the  image  bands.  Within  each 
band,  however,  a  narrowband  assumption  is  usually  made.  That  is,  the  band  center 
wavelength  defines  the  spectral  dependence.  The  parametric  modeling  tools,  however, 
explicitly  integrate  over  the  sensor  bandwidth. 

Most  modules  are  also  based  on  a  linear,  space  invariant  assumption.  For  example, 
many  effects  are  characterized  by  an  optical  transfer  function  (OTF)  which  operates  on  the 
spatial  Fourier  transform  of  an  image.  Nonlinear  and  space  variant  effects,  however,  are 
not  necessarily  incompatible  with  the  basic  model  formulation. 


3 


INPUT  IMAGE  EXPOSURE  DATA 

CONVERSION  GENERATION 


Figure  2-1:  Image  Simulation  Mode 


IMAGING  SPECTRAL  DATA  PERFORMANCE  l  PERFORMANCE 
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Figure  2-2:  Parametric  Performance  Mode 


Table  2-1:  IBSM  Modules 


ATMOSPHERIC  EFFECTS 

PROCESSING  EFFECTS 

MODTRAN  Input 

Image  Interlace 

Path  Trans/Rad 

Turbulence  OTF 

Aggregation 

Aero-Optic  OTF 

DATA  LINK  EFFECTS 

SENSOR  OPTICS  EFFECTS 

JPEG  Compress 

JPEG  Decompress 

Diffraction  OTF 

TSVQ  Compress/Decompress 

Defocus  OTF 

Wavelet  Comp/Decomp 

Jitter  OTF 

Finite  Buffer 

Drift  OTF 

Wavefront  OTF 

Channel  Errors 

User  OTF  (1-D) 

User  OTF  (2-D) 

PARAMETRIC  PERFORMANCE 

Radiometry 

Create  Spectra 

Create  MTF 

E/O  DETECTOR  EFFECTS 

Sensor  Performance 

Exposure  Control 

Detector  Size  OTF 

TDIOTF 

Image  Quality 

CTEOTF 

Diffusion  OTF 

UTILITIES 

Detection 

VIFFtoIBSM 

Noise 

List  IBSM  File  Data 

Quantization 

Transform 

Sampling 

Create  Plot  Data 

Nonuniformity 

Create  Radiance  Image 
Ground/Angle  Transform 

Load  Variables 

XGraph  Plot 
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2.2  MODELING  ENVIRONMENT 


IBSM  has  been  upgraded  from  an  earlier  version  [1]  to  operate  as  a  set  of  toolboxes 
within  the  Khoros  2.x  Cantata  environment.  Cantata  is  an  image  and  signal  processing 
prototyping  environment  which  is  widely  used  throughout  the  government  and  university 
research  communities.  What  Cantata  provides  is  the  flexible  modeling  environment  and 
graphical  user  interface  (GUI)  which  makes  IBSM  relatively  easy  to  use.  In  addition. 
Cantata  provides  a  large  set  of  standard  image  and  signal  processing  and  display  utilities 
which  are  at  the  user's  disposal. 

The  IBSM  software  consists  of  several  functional  layers  of  code  written  in  the  C 
programming  language.  The  executive  layer  was  written  using  the  Khoros  development 
system,  and  includes  all  the  Khoros-specific  functions,  including  file  I/O  and  graphical  user 
interface.  The  Khoros  library  function  layer  consists  of  the  interface  between  the  Khoros 
executive  code  and  the  sensor  modeling  library,  which  contains  the  fundamental  modeling 
routines.  The  sensor  modeling  library  is  a  set  of  C  subroutines  which  has  been  written 
independent  of  Khoros  or  any  other  executive  code.  Numerical  Recipes  and  MODTRAN 
are  used  for  analytical  computations  and  atmospheric  modeling. 

2.3  CAPABILITIES  SUMMARY 

IBSM  currently  consists  of  six  toolboxes:  Atmosphere,  Optics,  Detector,  Data 
Link,  Performance,  and  Utilities.  Each  of  these  toolboxes  are  briefly  described  here.  A 
more  detailed  description  is  given  in  Section  3. 

The  Atmosphere  toolbox  contains  the  tools  for  simulating  image  propagation  along 
the  atmospheric  path  between  the  imaged  area  and  the  sensor,  including  atmospheric  path 
transmission  and  radiance,  long  path  turbulence,  and  aero-optical  boundary  layer 
turbulence.  Atmospheric  path  transmission  and  radiance  effects  are  computed  by  executing 
MODTRAN.  A  separate  module  was  written  to  provide  a  graphical  user  interface  to 
MODTRAN  and  allow  simplified  entry  of  all  the  defining  atmospheric  parameters. 

The  Optics  toolbox  contains  the  tools  for  characterizing  the  image  effects  of  an 
optical  sensor  system.  These  effects  are  primarily  modeled  as  optical  transfer  functions 
(OTFs)  and,  therefore,  generally  operate  on  a  spatial  spectrum  (Fourier  transform  of  an 
image)  rather  than  the  image  itself.  OTFs  are  included  for  aperture  diffraction,  defocus. 
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wavefront  irregularity,  line-of-sight  drift,  and  line-of-sight  jitter.  Modules  for 
incorporating  user  defined  OTF  data  in  the  form  of  one-dimensional  and  two-dimensional 
OTFs  or  point  spread  functions  (PSF)  are  also  available.  Finally,  a  radiometric  tool  is 
provided  for  transforming  a  pupil  plane  radiance  image  into  a  focal  plane  irradiance  image, 
including  stray  radiance  sources. 

The  Detector  toolbox  contains  the  tools  for  characterizing  the  electro-optical 
detection  process  at  the  focal  plane.  Once  again,  several  of  the  effects  are  modeled  as 
OTFs,  including  that  due  to  the  finite  detector  element  size,  the  time-delay-integrate  (TDI) 
process,  charge  transfer  inefficiency,  and  carrier  diffusion.  Conversion  of  the  focal  plane 
irradiance  image  to  a  sampled  electronic  image  is  performed  by  the  combination  of  a 
photoelectronic  detection  and  sampling  tool.  The  non-ideal  performance  of  the  detector  can 
be  characterized  by  noise,  quantization,  and  nonuniformity  modules. 

The  Performance  toolbox  contains  the  modules  written  primarily  for  parametric 
sensor  modeling.  This  includes  modules  for  generating  atmospheric  and  sensor  spectral 
profiles;  computing  MTF  from  an  OTF  chain;  computing  SNR,  GRD,  and/or  NIIRS 
directly  from  sensor  parameters;  computing  line  response  functions,  SNR,  and  NIIRS 
from  an  image;  and  configuring  a  sensor  for  proper  exposure  control. 

The  Processing  toolbox  contains  modules  not  provided  by  the  standard  Khoros 
library  for  image  processing.  This  specifically  includes  functions  for  combining  image 
interlaces  and  aggregating  pixels  and/or  bands. 

The  Data  Link  toolbox  contains  the  tools  for  introducing  the  effects  of  lossy  data 
compression,  data  underflows  and  overflows  due  to  finite  transmission  buffers,  and 
channel  bit  errors  associated  with  transmitting  sensor  image  data  across  a  real 
communication  link. 

The  Utilities  toolbox  contains  routines  for  image  processing,  file  format 
conversion,  and  also  for  transforming  images  into  spatial  spectra,  and  vice  versa. 

A  complete  list  of  IBSM  modules  is  given  as  Table  2-1. 


8 


3.0  DETAILED  DESCRIPTION 


This  section  contains  a  detailed  description  of  the  physical  basis  of  the  IBSM 
toolboxes.  Before  proceeding,  however,  some  common  nomenclature  will  be  defined. 

For  routines  which  operate  on  radiance  imagery,  Lin(x,y)  and  Lout(x,y)  will  represent  the 
image  radiance  as  a  function  of  spatial  coordinates  x  and  y.  Analogous  notation  will  be 
used  for  reflectance  imagery  p(x,y),  temperature  imagery  T(x,y),  irradiance  imagery 

E(x,y),  and  imagery  in  detected  photoelectron  units  N(x,y).  The  spatial  coordinates  are 
generally  specified  in  angular  (radian)  units  unless  otherwise  noted. 

Optical  transfer  functions  are  characterized  by  H(u,v)  and  operate  on  spatial  spectra 
Gin(u,v)  to  form  Gout(u,v)  such  that 

G„„,(n,v)  =//(«,  v)G,,(m,v)  (3-1) 

where  u,v  are  spatial  frequency  (rad'^)  dimensions.  In  general,  the  OTF  and  the  spatial 
spectra  are  complex  (real  and  imaginary)  two-dimensional  distributions.  The  units  of 
G(u,v)  are  generally  arbitrary  and  H(u,v)  is  dimensionless. 

In  the  case  of  multiband  images  or  spatial  spectra,  the  basic  functions  outlined  in 
this  section  operate  separately  on  each  band.  This  is  implied  in  the  descriptions  to  follow, 
particularly  where  the  operations  are  wavelength  dependent.  Each  image  band  is 
characterized  by  its  center  wavelength  Xo  and  spectral  bandwidth  AX. 

In  the  Khoros  implementation  of  IBSM,  image  data  is  passed  between  modules  via 
files.  Generally,  the  imagery  is  in  floating-point  format  and  scaled  to  certain  physical  units. 
In  addition,  ancillary  information  is  contained  in  the  files  such  as  scaling  information, 
spectral  band  centers  and  bandwidths,  and  data  type  information.  This  is  included  in  the 
Khoros  file  structure  as  attributes.  The  need  for  this  additional  information  in  the  input  file 
is  denoted  for  each  of  the  modules  described. 
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3.1  ATMOSPHERE 


3.1.1  MODTRAN  Input 

This  function  provides  a  user  interface  for  generating  a  MODTRAN  path  data  file  as 
an  input  to  IBSM  functions  that  execute  MODTRAN,  including  Path  Trans/Rad,  Create 
Spectra,  and  Create  Radiance  Image.  The  output  file  is  in  ASCII  format,  as  given  in  Table 
3-1 .  All  MODTRAN-related  functions  have  been  upgraded  for  compatibility  with 
MODTRAN3. 

The  user  is  referred  to  the  applicable  MODTRAN  documents  [2,  3]  for  an  overview 
of  the  model  inputs.  A  few  comments  regarding  commonly  used  defaults  and  overrides  are 
given  here.  The  visibility  is  determined  from  the  haze  model  if  a  zero  value  is  entered, 
otherwise  it  is  overridden  by  the  user  entry.  The  surface  albedo  can  be  explicitly  entered 
(assumed  spectrally  flat),  or  determined  from  one  of  the  spectral  profiles  which  exist. 
Finally,  the  target  zenith  angle  is  really  measured  from  zenith  (MODTRAN  definition)  such 
that  -180°  is  a  downlooking  sensor. 

AU  modules  which  are  based  on  MODTRAN  require  a  path  data  file  and  an  optional 
working  directory  as  an  input.  The  module  creates  a  temporary  working  directory  (if  none 
is  specified);  copies  the  DIRAC,  sunl,  sun2,  and  refbkg  files  to  the  working  directory; 
forms  the  TAPES  input  file;  executes  MODTRAN3  (location  defined  by 
MODTRAN_HOME  environment  variable);  and  extracts  the  desired  data  from  the 
generated  TAPE6  or  TAPE?  files. 

3.1.2  Path  Trans/Rad 

This  function  applies  the  effects  of  atmospheric  path  transmission  and  radiance 
globally  to  an  input  radiance  image.  This  is  performed  in  two  steps.  First,  the  atmospheric 
transmission  x(X,)  and  path  radiance  Lpath  (A,)  as  a  function  of  wavelength  are  computed  by 
MODTRAN  based  on  the  path  data  file  created  by  the  MODTRAN  Input  function.  Then, 
the  band-average  transmission  and  path  radiance  is  applied  to  the  image 

1  *  \  P-2) 

AA  AA 
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Table  3- 1 :  Format  of  MODTRAN  Path  Data  File 


MODEL 

integer 

Model  atmosphere  (1-6) 

MULT 

integer 

Multiple  scattering  (0  =  off,  1  =  on) 

TBOUND 

float 

Surface  temperature  (K,  0.0  =  first  layer  temp.) 

SALE 

float 

Surface  albedo  (0-1,  -1,  -2,  -3,  -4,  -5,  -6) 

IHAZE 

integer 

Aerosol  haze  model  (0-6,  8-10) 

ICLD 

integer 

Cloud  model  (0-10, 18-20) 

VIS 

float 

Visibility  (km,  0.0  =  default) 

RAINRT 

float 

Rain  rate  (mm/hr) 

GNDALT 

float 

Ground  altitude  (km) 

IDAY 

integer 

Day  of  the  year  (1-365) 

ISOURC 

integer 

Celestial  source  (0  =  sun,  1  =  moon) 

ANGLEM 

float 

Moon  phase  angle  (deg) 

HI 

float 

Target  altitude  (km) 

H2 

float 

Sensor  altitude  (km) 

ANGLE 

float 

Sensor  zenith  angle  from  target  (deg) 

PARMl 

float 

Sun  azimuth  angle  from  sensor  LOS  (deg) 

PARM2 

float 

Sun  zenith  angle  from  target  (deg) 

IVl 

integer 

Minimum  wavenumber  (cm*l) 

IV2 

integer 

Maximum  wavenumber  (cm'l) 

IDV 

integer 

Wavenumber  sampling  (cm'l) 

IFWHM 

integer 

Wavenumber  resolution  (cm‘l) 
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It  is  imperative  that  the  wavelength  limits  specified  in  the  MODTRAN  path  data  file  cover 
each  spectral  band  specified  for  the  input  image  (attribute).  Spatial  scaling  attributes  for  the 
input  image  are  ignored.  The  input  image  must  be  in  radiance  (W/m^sr)  units. 

3.1.3  Turbulence OTF 

This  function  applies  the  effects  of  long  path  atmospheric  turbulence  to  an  image.  It 
is  characterized  as  an  OTF  and,  therefore,  operates  in  the  spatial  spectrum  domain.  The 
OTF  is  given  by  [4] 


H{u  v)  = 


(3-3) 


where  D  is  the  effective  aperture  diameter, 


p  =  ^u^+v^  (3-4) 

and  ro  is  the  correlation  diameter  defined  by  Fried  [5,  6]. 


The  correlation  diameter  can  be  either  user  specified  or  computed  from  an  index 
structure  parameter  model  along  the  line-of-sight  path.  In  the  latter  case  [7] 


To  =2.1 


5.%Ak 


Aq  cos  0 


..  (« 


K-K 


-1-3/5 


dh 


(3-5) 


where  hi  is  the  target  altitude,  hi  is  the  sensor  altitude,  ({)  is  the  target  nadir  angle,  and 
Cn^(h)  is  the  vertical  index  structure  parameter  distribution. 

Three  Cn^(h)  models  are  implemented,  all  of  which  are  parametrized  by  the  index 
structure  parameter  at  h  =  1  meter.  This  generally  ranges  from  lO  ^^  m'2/3  representing 
fairly  turbulent  conditions  to  lO'^'^m'^^  representing  fairly  benign  conditions  [8].  The  first 
model  is  a  simple  model  [7] 


C/(/i  =  lm) 
C„\h  =  \m) 


h 


h<lm 

h>lm 


(3-6) 
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The  second  is  a  model  proposed  by  Hufnagel  [9]  for  medium  to  high  altitudes,  and  is 
extended  to  lower  altitudes  using  a  power  coefficient  to  match  the  specified  Cn^(h)  at  h  =  1 
m[7].  The  functional  form  is 


C„\h)  = 


C^{h  =  \m)h  '■ 


C„2(/z  =  lm) 

15.85+log  (h-lm) 
3 


+  2.7xl0-''’e 


16^-fc/1500 


h<lm 


\m<h<  \km  (3-7) 
h>\km 


where  v  is  the  upper  atmospheric  wind  speed  (generally  in  the  18-36  m/s  range).  The  final 
model  is  an  approach  based  on  empirical  data  [10] 


C„\h  =  lm) 

h<lm 

C„\h  =  lm)h-^^ 

h  <  lOm 

h  <  1km 

3.16C/(/i  =  lm)r"' 

h  <  20km 

0 

h  >  20km 

(3-8) 


Figures  3-1  to  3-3  illustrate  the  index  structure  parameter  models  for  a  range  of 
conditions:  benign  (Cn^  =  lO'i^  m-2/3  at  h  =  1  m),  moderate  (Q^  =  lO'i^  m-2/3  at  h  =  1 
m),  and  turbulent  (Cn^  =  10'13  m'2/3  at  h  =  1  m).  For  the  Hufnagel  model,  upper 
atmosphere  wind  speeds  of  18, 27,  and  36  m/s  were  used,  respectively. 

The  a  parameter  in  Eq.  (3-3)  specifies  whether  the  OTF  is  of  the  long  exposure  (a 
=  0)  or  short  exposure  (a  =  0.5)  variety.  This  is  determined  by  the  temporal  nature  of  the 
turbulence  relative  to  the  sensor  integration  time.  In  general,  the  OTF  is  modeled  as  a  mix 
of  the  long  exposure  Hl  (u,v)  and  short  exposure  Hs  (u,v)  OTFs  given  by 

H(m,v)  =  +  /fi(M,v)(l  -  )  (3-9) 

where  Tj  is  the  pixel  dwell  (or  integration  time)  and  V  is  the  effective  velocity  of  the 
atmosphere  relative  to  the  sensor  pupil  plane  axes. 
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Hufnagel  Model 


Cn2  (m-2/3) 


Stable 

moderate 


le+00  le-K)l  le+02 


Altitude  (km) 


Figure  3-2:  Hufnagel  Index  Structure  Parameter  Model 
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Cn2(m-2/3) 


Empirical  Model 


Figure  3-3:  Empirical  Index  Structure  Parameter  Model 
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Spatial  scaling  (rad  ^  units)  and  spectral  band  attributes  must  be  available  for  the 
input  image,  and  the  data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced 
from  an  image  using  the  Transform  tool. 

3.1.4  Aero-Optic  OTF 

This  function  applies  the  effects  of  boundary  layer  aero-optical  turbulence  to  an 
image.  It  is  characterized  by  an  OTF  and,  therefore,  operates  in  the  spatial  spectrum 
domain.  The  OTF  is  given  by  [1 1] 


H(u,v)  =  e  L  J  (3.10) 

where  a(j,2  is  the  turbulence-induced  wavefront  phase  variance  and  Lx,  Ly  are  the  spatial 
correlation  lengths  of  the  turbulent  boundary  layer. 

The  phase  variance  is  given  by 

(3-11) 

Aq 

where  G  is  the  Gladstone-Dale  coefficient  (0.22  x  lO’^  m^/kg),  p'  is  the  rms  density 
fluctuation,  Lz  is  the  correlation  length  of  the  boundary  layer  in  the  normal  direction,  and 
tbi  is  the  boundary  layer  thickness. 

The  model  is  based  on  a  "rule-of-thumb"  equation  for  the  rms  density  fluctuation 
and  correlation  lengths  [12].  Specifically, 

p'  =  0.2[po-p]  (3-12) 


where. 


P  = 


_ Po _ 

l  +  r(7-l)MV2 


(3-13) 


17 


Po  is  the  free  stream  air  density,  r  =  0.89, 7=  1.4,  and  M  is  the  Mach  number  (v/305  m/s). 
The  correlation  lengths  are  approximated  by 


II 

p 

(3-14) 

II 

II 

p 

N> 

(3-15) 

The  boundary  layer  is  then  specified  by  its  thickness  tbi  along  with  the  free  stream  air 
density  and  platform  velocity.  Table  3-2  provides  a  typical  free  stream  air  density  as  a 
function  of  altitude. 


Spatial  scaling  (rad*^  units)  and  spectral  band  attributes  must  be  available  for  the 
input  image,  and  the  data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced 
from  an  image  using  the  Transform  tool. 

3.2  IBSM  OPTICS 

3.2.1  Diffraction  OTF 


This  function  applies  the  effects  of  aperture  diffraction  to  an  image.  It  is 
characterized  by  an  OTF  and,  therefore,  operates  in  the  spatial  spectrum  domain.  For  an 
unobscured  circular  aperture  [13] 


H{u,v)  = 


cos 


-1 


VPoJ 


Po 


T 

K.P0J 


P^Po 

P>Po 


(3-16) 


where  p  is  defined  by  Eq.  (3-4)  and 


(3-17) 


where  D  is  the  effective  aperture  diameter.  For  a  rectangular  aperture  [13] 
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Table  3-2;  Typical  Free  Stream  Air  Density  as  a 
Function  of  Attitude 


Altitude  (kmi 


tir  Uensit 
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H{u,v)  =  K 


\^xj 


ydyj 


(3-18) 


where  dx,  dy  are  the  aperture  widths  in  their  respective  directions  and 


A(0  = 


fi-lil  III  SI 
I  0  |^|>1 


(3-19) 


In  the  case  of  an  obscuration,  only  the  circular  symmetric  case  is  modeled 
(obscuration  input  is  ignored  for  rectangular  aperture).  In  this  case,  the  OTF  is  given  by 
[14] 


H(u,v)  = 


A  +  B  +  C 

i-v" 


(3-20) 


A= 


cos 


-All- 

Po 


T 

vPoy 


P^Po 

P>Po 


(3-21) 


B= 


27]^ 


n 


cos 


-1 


(_R_] 


^  p 


VPo 


1- 


mj 


P^VPo 

p>m 


(3-22) 
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277  .  , 

— -  -  Sin  <p  + 
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-Irf 
^1  +  11^^ 


\  K  j 


- ^tan 

n 

0 


vl-^y 


(j)\ 

tan  - 

2) 


1-7] 

P<^Po 


l-T]  ^  .1+t? 

-ir^Po^P^—^Po 


l+Tl 

P>-T-^Po 


(3-23) 


(j)  =  cos 


-1 


\  +  j]^-(2p/pof 

277 


(3-24) 
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p  is  defined  in  Eq.  (3-4),  po  =  D/Xo,  and  q  is  the  relative  obscuration  (ratio  of  obscuration 
to  aperture  diameters). 

Spatial  scaling  (rad'*  units)  and  spectral  band  attributes  must  be  available  for  the 
input  image,  and  the  data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced 
from  an  image  using  the  Transform  tool. 

3.2.2  Defocus OTF 

This  function  applies  the  effects  of  optics  defocus  or  blur  to  an  image.  It  is 
characterized  by  an  OTF  and,  therefore,  operates  in  the  spatial  spectrum  domain.  The 
defocus  (or  blur)  function  is  modeled  as  a  gaussian  distribution  for  which  the  1/e  blur 
radius  in  angular  space  is  specified  independently  for  the  x  and  y  axes.  Mathematically 
[15] 


H{u,v)  =  e  ^ 


(3-25) 


where  Wx  and  Wy  are  the  1/e  blur  spot  radii  in  x  and  y. 

Although  not  exact,  this  module  can  be  used  to  model  OTF  degradations  due  to 
axial  misfocus.  The  actual  blur  function  will  be  dependent  on  the  aperture  characteristics, 
but  a  gaussian  blur  function  is  a  fair  approximation.  For  an  axial  defocus  defined  in  the 
object  domain  (focus  set  to  Rq,  object  at  R) 


w^y  =  0.622)1 


R  R, 


For  an  axial  defocus  Az  defined  in  the  image  domain. 


0.62DAZ 


w  =  • 


r 


(3-26) 


(3-27) 


Spatial  scaling  (rad'*  units)  attributes  must  be  available  for  the  input  image,  and  the 
data  must  correspond  to  a  spatial  spectmm,  which  can  be  produced  from  an  image  using 
the  Transform  tool. 
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3.2.3  Jitter OTF 


This  function  applies  the  effects  of  wideband  (with  respect  to  the  sensor  integration 
time)  jitter  of  the  optical  sensor  line-of-sight  to  an  image.  It  is  characterized  by  an  OTF 
and,  therefore,  operates  in  the  spatial  spectrum  domain.  By  virtue  of  the  wideband 
assumption,  the  jitter  is  modeled  as  a  zero  mean,  gaussian  random  process  resulting  in  the 
gaussian  expected  value  for  the  OTF  [16]: 

=  (3-28) 

where  Gx  and  Gy  are  the  rms  jitter  amplitudes  in  their  corresponding  angular  directions. 

Spatial  scaling  (rad'i  units)  attributes  must  be  available  for  the  input  image,  and  the 
data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced  from  an  image  using 
the  Transform  tool. 

3.2.4  Drift  OTF 

This  function  applies  the  effects  of  constant  angular  drift  of  a  sensor  line-of-sight  to 
an  image.  It  is  characterized  by  an  OTF  and,  therefore,  operates  in  the  spatial  spectrum 
domain.  The  OTF  is  defined  by 


H(u,  v) = sinc(a^M)sinc(ay  v) 


(3-29) 


where 


=  (3-30 

and  ax  and  ay  are  the  angular  drift  magnitudes  over  the  sensor  integration  time  in  their 
corresponding  directions. 

Spatial  scaling  (rad-i  units)  attributes  must  be  available  for  the  input  image,  and  the 
data  must  correspond  to  a  spatial  spectmm,  which  can  be  produced  from  an  image  using 
the  Transform  tool. 
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3.2.5  Wavefront  OTF 


This  function  applies  the  effects  of  a  spatially  random  wavefront  error  to  an  image. 
It  is  characterized  by  an  OTF  and,  therefore,  operates  in  the  spatial  spectrum  domain.  The 
OTF  is  given  by  [17] 


H(u,v)  =  e  L  J  (3.31) 

where  the  wavefront  error  is  defined  in  the  sensor  pupil  plane  by  the  phase  variance  a((,2 
and  X  and  y  correlation  lengths  1^  and  Ly  (defined  in  spatial  aperture  dimensions).  The 
phase  variance  is  specified  by  an  rms  wavefront  error  (waves)  at  a  reference  wavelength 
and  scaled  to  the  image  band  center  wavelength(s)  in  inverse  proportion  to  the  wavelength. 

Spatial  scaling  (rad’i  units)  and  spectral  band  attributes  must  be  available  for  the 
input  image,  and  the  data  must  correspond  to  a  spatial  spectmm,  which  can  be  produced 
from  an  image  using  the  Transform  tool. 

3.2.6  User  OTF  (1-D) 

This  function  applies  an  OTF  specified  by  one-dimensional  user  data  to  an  image. 
As  it  is  an  OTF,  it  operates  in  the  spatial  spectrum  domain.  The  data  is  contained  in  an 
ASCII  file  in  formats  accommodating  real  and  complex  data  specified  in  radial  form 
(radial),  with  separability  in  the  x  and  y  axes  (x  and  y),  or  with  separability  and  similarity 
in  the  x  and  y  directions  (x  or  y).  The  user  data  is  interpolated  using  linear  or  spline 
techniques  to  the  input  spatial  frequency  grid  prior  to  forming  the  2-D  OTF. 

The  various  file  formats  characterizing  the  user  OTF  are  given  in  Table  3-3.  In  the 
case  of  radial  data,  the  OTF  is  interpolated  onto  the  (u,v)  grid  at  each  point  relative  to  its 
radius  defined  by  Eq.  (3-4).  For  x  and  y  data,  the  input  data  is  interpolated  onto  two 
functions  Hx(u)  and  Hy(v),  and  the  two-dimensional  OTF  is  defined  by 

H{u,v)  =  HMtiy(y)  (3-32) 
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Table  3-3:  User  OTF  (1-D)  File  Formats 


SF(1) 


SF(N) 


OTFreal(l)  [OTFimagd)] 

OTFreal(N)  [OTFimagCN)] 

(a)  Radial  Data 


SF(1) 


SF(N) 


OTFreal(l)  [OTFimag(l)] 


OTFreal(N)  [OTFiinag(N)] 

(b)  X  and  Y  Data 


SFx(l)  OTFx,real(l)  [OTFx,imag(l)]  SFyd)  OTFY,reald)  [OTFY,iinagd )] 


SFx(N)  OTFx,real(N)  [OTFx,imag(N)]  SFy(N)  OTFY,real(N)  [OTFY,imag(N)] 

(c)XorYData 


[]  =  data  element  exist  only  for  complex  data 
N  lines  of  space  delimited  data  in  all  cases 
SF  =  Spatial  Frequency  (rad'l) 

OTF  =  Optical  Transfer  Function 
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For  jc  or  y  data,  Hx(u)  and  Hy(v)  are  identical  by  definition.  Spatial  scaling  (rad  ^  units) 
attributes  must  be  available  for  the  input  image,  and  the  data  must  correspond  to  a  spatial 
spectmm,  which  can  be  produced  from  an  image  using  the  Transform  tool. 

3.2.7  UserOTF(2-D) 

This  function  applies  an  OTF  specified  by  two-dimensional  user  data  to  an  image. 
As  it  is  an  OTF,  it  operates  in  the  spatial  spectmm  domain.  The  user  data  containing 
scaling  information  is  contained  in  an  image  file.  This  data  can  correspond  to  an  OTF 
(rad'i  scaling  attributes)  or  a  point  spread  function  (PSF,  rad  scaling  attributes). 

In  the  case  of  OTF  data,  the  user  defined  data  is  interpolated  onto  the  input  spatial 
frequency  grid  using  a  bilinear  technique.  In  the  case  of  PSF  data,  the  user  OTF  is 
computed  at  each  grid  point  by 


H{u,v)  =  —  E  E  /z  (i,j)  (3-33) 

Ik 

where  Xi  and  yj  are  the  angular  positions  corresponding  to  the  PSF  indices.  If  specified, 
the  OTF  data  is  normalized  to  unity  at  zero  spatial  frequency  in  either  case. 

Spatial  scaling  (rad'*  units)  attributes  must  be  available  for  the  input  image,  and  the 
data  must  correspond  to  a  spatial  spectmm,  which  can  be  produced  using  the  Transform 
tool. 


3.2.8  Radiometiy 


This  function  performs  a  radiometric  conversion  of  an  effective  pupil  plane  radiance 
image  to  a  focal  plane  irradiance  image,  including  the  sensor  throughput,  thermal  radiation 
of  the  optical  chain,  and  stray  radiance: 


['^opAx,y)+  e,p,LBs{X,T^^,)AX  +  L, 


(3-34) 


where  (f/#)  is  the  effective  sensor  F-number,  x  is  the  optical  train  transmission  including 
obscuration,  Eopt  is  the  optical  train  effective  emissivity,  Topt  is  the  optical  train  effective 
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temperature,  Ls  is  the  in-band  stray  radiance  (non-thermal),  and  Lbb  is  a  blackbody  spectral 
radiance  function 


=  (3-35) 

Spectral  scaling  attributes  must  be  available  in  the  input  file  header,  and  the  image 
data  must  be  in  radiance  (W/m^sr)  units.  The  output  image  data  is  in  irradiance  (W/m^) 
units. 


3.3  IBSM  DETECTOR 


3.3.1  Detector  Size  OTF 


This  function  applies  the  effects  of  the  finite  detector  width  to  an  image.  It  is 
characterized  by  an  OTF  and,  therefore,  operates  in  the  spatial  spectrum  domain.  The 
detector  is  assumed  to  have  rectangular  form,  which  generates  an  OTF  with  a  sine 
distribution  [18].  It  is  extended,  however,  to  include  pixel  aggregation  effects: 


H(u,v)  =  sine 


K  f  J 


sme 


V  /  y 


N^-l  N,-l 

E  E  cos  (bi ; 

i=0  j=0  '’■> 


(3-36) 


where  w^  and  Wy  are  the  detector  sizes  in  x  and  y,  f  is  the  optical  system  focal  length,  Nx 
and  Ny  are  the  number  of  aggregated  pixels  in  x  and  y. 


<l).j  =  Ik 
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(3-37) 


and  Px  and  Py  are  the  detector  center-to-center  spacings  (pitch)  in  x  and  y.  The  summation 
in  Eq.  (3-36)  results  from  the  sum  of  aggregated  detector  elements,  with  a  phase  shift 
corresponding  to  the  relative  spacing  between  elements  in  the  image  domain. 

Spatial  scaling  (rad'i  units)  attributes  must  be  available  for  the  input  image,  and  the 
data  must  correspond  to  a  spatial  spectmm,  which  can  be  produced  from  an  image  using 
the  Transform  tool. 
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3.3.2  TDIOTF 


This  function  applies  the  effects  of  time-delay-integration  (TDI)  to  an  image.  It  is 
characterized  by  an  OTF  and,  therefore,  operates  in  the  spatial  spectrum  domain.  The  TDI 
can  occur  in  either  the  x  or  y  directions,  and  may  exhibit  a  mismatch  between  the  TDI 
clocking  rate  and  the  image  motion  rate.  The  OTF  is  given  by: 


H{u,v)  =  sine 
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(3-38) 


where  w  is  the  spatial  frequency  in  the  TDI  direction  (u  or  v),  d  is  the  detector  width  in  the 
TDI  direction,  f  is  the  optical  system  focal  length,  Ntdi  is  the  number  of  TDI  stages, 

Nphases  is  the  number  of  clock  phases  per  transfer,  and  P  is  the  ratio  of  the  TDI  clocking 
rate  to  image  motion  rate.  The  sine  function  results  due  to  the  image  motion  over  each  TDI 
stage,  while  the  summation  results  from  the  relative  shifts  between  TDI  stage  images  which 
arises  due  to  the  rate  mismatch. 


Spatial  scaling  (rad'i  units)  attributes  must  be  available  for  the  input  image,  and  the 
data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced  from  an  image  using 
the  Transform  tool. 

3.3.3  CTEOTF 


This  function  applies  the  effects  of  charge  transfer  inefficiency  of  a  charge  coupled 
device  (CCD)  detector  array  readout  to  an  image.  It  is  characterized  by  an  OTF  and, 
therefore,  operates  in  the  spatial  spectrum  domain.  The  OTF  is  given  by  [19] 

H{u  v)  =  P9 

where  Nx  and  Ny  are  the  number  of  transfers  in  x  and  y,  Px  and  py  are  the  center-to-center 
detector  spacings  (pitch)  in  x  and  y,  f  is  the  optical  system  focal  length,  Nphase  is  the 
number  of  clock  phases  per  transfer,  and  e  is  the  charge  transfer  efficiency  (per  phase). 

Spatial  scaling  (rad-i  units)  attributes  must  be  available  for  the  input  image,  and  the 
data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced  from  an  image  using 
the  Transform  tool. 
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3.3.4  Diffusion  OTF 


This  function  applies  the  effects  of  minority  carrier  diffusion  in  a  charge  coupled 
device  (CCD)  sensor  to  an  image.  It  is  characterized  by  an  OTF  and,  therefore,  operates  in 
the  spatial  spectrum  domain.  The  OTF  is  given  by  [20] 


H{u,v) 


1  +  aAL[p) 

^-aALD 

1  +  cxALO 


(3-40) 


where 


AL{p)  = 
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(3-41) 


p  is  given  by  Eq.  (3-4),  a  is  the  carrier  spectral  diffusion  coefficient,  ALD  is  the  depletion 
layer  width,  ALO  is  the  diffusion  length,  and  f  is  the  optical  system  focal  length.  Table  3-4 
provides  the  carrier  spectral  diffusion  coefficient  for  silicon  at  various  wavelengths. 

Spatial  scaling  (rad’^  units)  attributes  must  be  available  for  the  input  image,  and  the 
data  must  correspond  to  a  spatial  spectmm,  which  can  be  produced  from  an  image  using 
the  Transform  tool. 

3.3.5  Detection 


This  function  performs  a  photoelectronic  conversion  of  a  focal  plane  irradiance 
image  to  a  detected  photoelectron  image,  including  detector  sensitivity,  aggregation,  dark 
current,  and  saturation.  The  conversion  process  is  given  by 


E{x,y)  -h - 

he  3^  3^  g 


(3-42) 


where  r\  is  the  quantum  efficiency,  Tj  is  the  pixel  dwell  time,  Ntdi  is  the  number  of  TDI 
stages,  Nx  and  Ny  are  the  number  of  aggregated  pixels  in  x  and  y,  Wx  and  Wy  are  the 
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Table  3-4:  Carrier  Spectral  Diffusion  Coefficient  for 
Silicon  at  Various  Wavelengths 
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detector  widths  in  x  and  y,  Jd  is  the  dark  current  density,  h  is  Planck's  constant,  c  is  the 
speed  of  light  in  a  vacuum,  and  q  is  the  electronic  charge. 

Spectral  scaling  attributes  must  be  available  for  the  input  image,  and  the  image  data 
must  be  in  irradiance  (W/m^  units).  The  output  image  data  is  in  units  of  photoelectron 
counts. 


3.3.6  Noise 


This  function  adds  a  random  noise  process  to  an  input  image.  Two  types  of  noise 
are  modeled.  The  first  is  signal  dependent,  Poisson  distributed  shot  noise,  which  can  be 
either  enabled  or  disabled.  The  second  is  signal  independent,  additive  noise  that  is 
assumed  zero  mean,  gaussian  distributed  with  a  specified  variance.  The  noise  sources  are 
assumed  both  jointly  and  pixel-to-pixel  independent.  The  input  and  output  image  data  must 
be  in  photoelectron  units. 


3.3.7  Quantization 

This  function  quantizes  a  detected  photoelectron  image  through  a  scalar  roundoff 
process.  Mathematically, 
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(3-43) 


Ku,{x^y)  = 
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where  N^in  and  Nmax  are  the  minimum  and  maximum  input  levels  of  the  quantizer  and  b  is 
the  number  of  digitization  levels  (bits). 

The  input  image  data  scaling  is  arbitrary,  but  must  be  consistent  with  the  minimum 
and  maximum  levels. 


3.3.8  Sampling 

This  function  samples  an  input  image  on  a  rectilinear  grid  with  arbitrary  spacing. 
The  sample  locations  in  angular  (radian)  units  are  given  by 
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(3-44) 


^  =  ^min+('  +  ^Jy 

y  =  >’min+(y  +  ^,)y  (3-45) 

where  px  and  py  are  the  detector  center-to-center  spacings  in  x  and  y,  f  is  the  optical  system 
focal  length,  6x  and  5y  are  the  user  defined  relative  offset  of  the  sampling  grid  in  x  and  y 
(relative  to  pitch),  x^in  and  y^in  are  the  minimum  field  angles  of  the  input  image,  and  i  and 
j  are  indices  which  range  from  zero  to  the  maximum  values  such  that  the  angular  sample 
locations  do  not  exceed  the  input  image  maximum  field  angles  x^ax  and  ymax-  The 
resampling  process  is  achieved  using  a  bilinear  interpolation  of  the  input  image  data. 

Spatial  scaling  (rad  units)  attributes  must  be  available  for  the  input  image.  The 
output  image  scaling  attributes  will  vary  depending  on  the  sampling  parameters,  but  will 
likely  be  different  than  the  input. 

3.3.9  Nonuniformity 

This  function  applies  the  effects  of  random  detector  gain  and/or  offset 
nonuniformity  to  an  input  image.  The  nonuniformity  can  be  characterized  as  an 
independent  random  variable  between  image  samples,  or  the  spatial  correlation 
characteristics  can  be  specified  by  an  input  autocorrelation  image  file.  In  both  cases, 
realizations  of  gaussian  distributed  random  gain  and  offset  patterns  are  generated 
(independent  random  sample  for  each  input  image  pixel).  The  gain  pattern  G(x,y)  is  unity 
mean  with  a  standard  deviation  given  by  the  gain  nonuniformity  and  the  offset  pattern 
0(x,y)  is  zero  mean  with  a  standard  deviation  given  by  the  offset  nonuniformity  times  the 
maximum  signal  level. 

When  a  spatial  autocorrelation  file  is  given,  it  is  first  resampled  into  a  convolution 
kernel  with  an  x  and  y  sample  spacing  identical  to  the  input  image.  This  is  performed  by  a 
bilinear  interpolation.  Next,  it  is  normalized  to  exhibit  unit  variance.  Finally,  the  gain  and 
offset  patterns  are  convolved  by  it  to  generate  the  desired  spatial  correlation  characteristics. 

Once  the  appropriate  random  gain  and  offset  patterns  are  generated,  the  output 
image  is  computed  by 
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L,{x^y)  =  G{x,y%{x,y)  +  0{x,y) 


(3-46) 


The  units  of  the  input  image  are  arbitrary,  but  must  correspond  to  the  specified  maximum 
signal  level  for  proper  scaling  of  the  offset  nonuniformity.  Spatial  scaling  attributes  must 
exist  for  the  input  and  autocorrelation  files  when  an  autocorrelation  file  is  used.  The  units 
are  arbitrary,  but  must  be  consistent. 

3.4  DATA  LINK 

3.4.1  JPEG  Compress 

This  function  compresses  an  8  or  12  bit  image  using  the  Joint  Photographic  Experts 
Group  compression  scheme  [21]  according  to  MIL-STD-188-198A  [22].  The  routine 
inputs  a  multiband  IBSM  image,  an  image  quality  factor,  and  a  restart  interval  length.  It 
compresses  each  band  separately  and  returns  both  a  compressed  data  stream  and  a 
multiband  bits/block  image.  Compression  is  performed  using  code  from  the  National 
Image  Transmission  Format  Standard  (NIFTS)  JPEG  Compression  (NJC)  Software 
Library  (NJC)  [23]. 

Bits/Block  Image  CBB  Image) 

The  BB  image  for  some  arbitrary  image  has  one  element  for  each  MCU  in  the 
image.  This  element  contains  the  number  of  bits  in  the  compressed  representation  of  that 
MCU,  and  is  used  as  an  input  for  the  Finite  Buffer  routine. 

The  BB  image  is  derived  by  counting  the  number  of  bits  in  each  restart  interval  and 
averaging  this  over  the  number  of  MCU’ s  in  a  restart  interval.  The  BB  image  is  set  equal 
to  this  average  for  the  entire  restart  interval.  Therefore  the  shorter  the  restart  interval 
length,  the  more  accurate  the  BB  image.  And  (nearly)  perfectly  accurate  BB  images  are 
produced  if  and  only  if  the  restart  interval  is  one  MCU  long. 
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Minimum  Coded  Unit  rMCU) 


MCU  is  a  JPEG  term  which  stands  for  Minimum  Coded  Unit.  To  define  it,  start 
with  any  single-band  image  and  divide  it  into  8x8  blocks  which  cover  the  whole  image  but 
yet  do  not  overlap.  Then  each  of  these  blocks  is  an  MCU. 

Image  Quality  Factor  (lOF) 

IQF  is  an  NIC  and  IBSM  term  which  denotes  a  number  from  1  to  5  that  controls 
the  accuracy  with  which  the  image  is  reconstructed.  Lower  values  mean  that  fewer  bits  are 
used  in  the  compressed  file.  Unfortunately  they  also  mean  that  the  reconstructed  image  is 
cruder.  This  range  of  values  (i.e.  1  to  5)  was  set  by  MIL-STD-188-198A,  but  is  not  part 
of  the  JPEG  standard  itself. 

Restart  Interval 


This  is  a  JPEG  term  which  denotes  a  collection  of  one  or  more  consecutive  MCU’s. 
After  the  compressed  data  for  each  restart  interval,  the  compressor  transmits  a  code.  This 
may  allow  the  receiver  to  regain  synchronization  when  a  bit  error  occurred  in  the 
compressed  data  for  that  restart  interval.  In  other  words,  channel  errors  should  be  confined 
to  the  restart  interval  in  which  they  occurred.  So  for  optimal  error  protection  we  want  the 
restart  intervals  as  short  as  possible.  The  tradeoff  is  that  each  restart  code  takes  about  20 
bits  to  transmit,  so  the  shorter  the  restart  interval  the  more  bits  are  spent  on  overhead. 

To  comply  with  MIL-STD-188-198A  an  image’s  restart  interval  must  be  no  larger 
than  the  number  of  MCU’s  it  takes  to  span  the  width  of  that  image. 

3.4.2  JPEG  Decompress 

This  function  is  the  inverse  of  JPEG  Compress.  It  inputs  a  JPEG  compressed  data 
stream  and  returns  an  approximation  to  the  original  image. 

3.4.3  TSVQ  Compress/Decompress 

Vector  quantizers  are  often  used  to  compress,  with  significant  loss,  sampled  data 
such  as  images  or  sound.  The  class  of  vector  quantizers  has  a  particular  subset  whose 
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members,  called  Tree-Stmctuied  Vector  Quantizers  (TSVQ’s),  can  achieve  a  given  level  of 
compression  performance  at  a  low  complexity.  [24]  describes  them  in  detail. 

This  function  has  two  modes.  The  first  reads  an  input  image  and  a  TSVQ 
codebook.  It  compresses  the  input  using  that  codebook  and  mainly  produces  a  sequence  of 
indices  that  represents  the  input  image.  The  second  mode  reads  a  sequence  of  indices  and 
outputs  the  decompressed  image. 

The  vector  quantizer  (VQ)  described  here  is  optimized  to  minimize  the  mean- 
squared  error  (MSE)  between  the  reconstmcted  image  and  the  original  image.  This  is  the 
customary  way  that  Vqs  are  designed,  mostly  because  it  is  conceptually  the  simplest 
distortion  measure.  Unfortunately  it’s  bad  to  use  MSE  distortion  if  the  reconstructed 
images  will  be  primarily  interpreted  by  humans,  although  this  may  work  ok  for  machine 
interpretation.  When  human  viewing  is  important,  VQ’s  (of  this  type)  are  usually  used  as 
components  in  fancier  data  compressors  such  as  wavelet  coders. 

This  TSVQ  program  is  based  on  programs  from  [25].  It  is  basically  a 
straightforward  TSVQ  encoder  and  decoder,  with  the  addition  of  restart  markers.  The 
restart  marker  scheme  is  similar  to  the  one  in  the  JPEG  standard  [21].  It  starts  by  picking 
some  natural  number  restart  interval  length  (RI_len),  then  takes  the  stream  of  (possibly 
variable-length)  indices  and  puts  a  marker  after  every  RI_len  of  them.  A  marker  consists  of 
two  bytes.  The  first  is  always  7F.  The  second  is  a  number  that  cycles  through  the 
markers. 

At  the  receiver,  the  data  is  parsed  into  segments  delimited  by  markers.  Each 
segment’s  marker  determines  the  location  in  the  image  to  begin  placing  the  pixels 
corresponding  to  that  segment’s  indices.  This  works  perfectly  if  the  markers  never  get 
corrupted.  One  complication  is  that  7F’s,  which  indicate  markers,  can  also  appear  naturally 
in  the  bits  composing  the  indices.  This  is  bypassed  by  stuffing  a  00  byte  after  each  natural 
occurrence  of  7F.  The  receiver  then  needs  to  remember  that  7F  00  needs  to  be  translated  to 
just  7F  before  the  decompressor  sees  it. 

Unfortunately,  the  markers  can  get  corrupted.  One  coping  strategy  is  to  remember 
the  index  (i.e.  the  position  in  markers  [])  of  the  last  decoded  marker.  Then  look  for  the 
next  marker  whose  index  is  greater.  If  its  index  isn’t  exactly  one  greater,  then  the  receiver 
knows  that  an  error  has  been  made  but,  for  simplicity,  doesn’t  try  to  estimate  the  missing 
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marker’s  location.  With  this  strategy,  if  an  error  occurs  in  a  marker,  then  the  entire 
segment  corresponding  to  that  marker  is  lost.  Also,  if  an  error  occurs  in  the  indices,  then  at 
most  the  rest  of  the  segment  is  lost. 

Bits/Unit  Time  File 


Imagine  a  mesh  whose  interstices  are  8x8  pixels  square.  Overlay  this  on  the  input 
image.  The  bits/(unit  time)  file  has  one  32-bit  integer  per  interstice.  Each  of  these  integers 
is  the  number  of  bits  emitted  per  unit  time  when  the  TSVQ  is  compressing  the  pixels  in  that 
interstice.  It  turns  out  that  TSVQ’s  emit  a  constant  number  of  bits  per  unit  time  so  this  file 
is  not  terribly  exciting,  but  it’s  still  sometimes  needed. 

3 .4.4  Wavelet  Compress/Decompress 

Wavelet  coding  [24, 26, 27, 30]  tries  to  approximate  the  compression  performance 
of  very  high-dimensional  vector  quantization  at  a  reasonable  complexity.  Suppose  one 
wants  to  compress  images  that  are  N  pixels  long.  The  points  in  N-dimensional  space  that 
represent  these  images  will  typically  be  concentrated  in  a  cloud  that  occupies  only  a  small 
fraction  of  the  volume  of  this  space.  The  basic  idea  of  wavelet  encoding  is  to  apply  an  N  x 
N  wavelet  transform  matrix  to  each  of  these  points.  This  effectively  orients  the  N 
coordinate  axes  with  the  cloud  so  that  it’s  easier  to  describe  the  location  of  a  particular 
image  point.  The  description  is  still  done  with  a  vector  quantizer  (VQ),  but  the  VQ  has 
much  less  work  because  it  happens  that  it  can  now  consider  many  fewer  than  N  dimensions 
at  a  time. 

The  wavelet  decoder  takes  the  indices  that  the  VQ  produced  and  reconstructs  the 
image  in  the  transformed  coordinate  system.  Finally,  it  returns  the  image  to  the  normal 
coordinate  system  by  multiplying  by  the  inverse  wavelet  transform.  For  a  simple  example 
of  the  advantages  of  a  good  coordinate  system,  consider  the  description  of  a  point’s 
location  on  a  particular  circle.  In  polar  coordinates  the  description  is  economical  because  it 
requires  just  one  number:  the  point’s  angle  with  respect  to  some  reference  ray. 

Conversely,  in  cartesian  coordinates  the  description  requires  two  numbers:  the  x  and  y 
positions. 

This  function  has  two  modes:  compress  and  decompress.  It  is  executed  based  on  a 
codebook  from  a  training  data  set  similar  to  that  with  TSVQ. 
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3.4.5  Finite  Buffer 


This  program  simulates  the  buffering  of  an  IBSM  format  image.  It  inputs  the  image 
and  its  corresponding  bits/block  image  (see  the  documentation  for  the  JPEG  glyphs  for 
more  information  about  these).  The  latter  is  produced  automatically  by  IBSM  compressor 
glyphs,  or  it  may  be  produced  manually  for  non-compressed  images.  The  outputs  of  this 
function  consist  of  the  buffered  data,  which  may  be  different  from  the  original  data  if 
overflows  and  underflows  occurred.  It  also  produces  a  graph  that  represents  how  many 
bits  were  in  the  buffer  at  each  time  step. 

This  program  keeps  track  of  two  kinds  of  events.  If  “Direction  Toggle”  is  0  then 
these  are  bit  arrivals  and  block  departures,  otherwise  they  are  bit  departures  and  block 
arrivals.  This  program  assumes  that  the  time  axis  has  been  partitioned  into  intervals  of 
length  “Step  size”.  It  examines  one  interval  at  a  time,  simulates  the  events  that  happen  in 
that  interval  and  then,  if  requested,  saves  the  number  of  bits  remaining  in  the  buffer  at  the 
end  of  the  interval. 

The  bookkeeping  it  does  is  conceptually  simple.  The  only  twist  is  that  there  is  no 
area  of  memory  that  serves  as  an  actual  buffer.  Instead,  the  program  keeps  track  of  the 
number  of  bits  in  an  imaginary  buffer.  If  an  incoming  bit  will  fit  in  this  buffer,  then  it  is 
placed  directly  on  the  output  since  if  there  really  were  a  buffer,  then  the  bit  would  get  there 
eventually  anyway.  This  is  simpler  than  first  writing  it  into  the  buffer,  and  then  writing  it 
into  the  output. 

The  buffer  model  is  based  on  a  fixed  rate  source  assumption;  that  is,  there  is  a  fixed 
number  of  seconds  between  arrivals  of  elements  (either  bits  or  blocks  dependent  on  the 
direction  toggle).  Both  buffer  overflows  and  underflows  can  be  simulated,  and  the  buffer 
can  be  initialized  to  an  initial  level  and/or  delay  prior  to  emptying. 

3.4.6  Channel  Errors 

This  program  inputs  an  IBSM  format  image  and  simulates  passing  it  through  a 
noisy  channel.  It  then  outputs  the  resulting  corrupted  image.  Any  combination  of  the 
following  channels  may  be  used:  binary  symmetric  (substitution)  channels  (BCS’s), 
insertion  channels,  deletion  channels,  and  burst  error  channels.  Each  bit  passing  through  a 
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BSC  has  the  same  probability  of  being  flipped.  An  insertion  channel  randomly  inserts  bits 
into  the  transmitted  data  stream,  while  a  deletion  channel  randomly  removes  bits.  Finally,  a 
burst  error  channel  corrupts  sequences  of  consecutive  bits. 

This  program  is  based  on  the  idea  that  for  each  bit  position  in  the  file  one  of  five 
events  can  happen:  (no  error,  bit  flip,  bit  insertion,  bit  deletion,  burst  error)  =  (e_0,  e_l, 
e_2,  e_3,  e_4).  Using  this  notation,  each  sequence  of  events  is  seen  to  be  an  element  of 
{e_0,  e_l,  e_2,  e_3,  e_4}*.  An  obvious  implementation  is  to  use  the  inverse 
transformation  method  (described  in,  for  example,  [30])  at  each  bit  position  to  map  the 
event  probabilities  into  a  sequence  from  this  set. 

This  program  uses  a  faster  alternative  that  is  based  on  the  equality  of  {  e_0,  e_l, 
e_2,  e_3,  e_4}*  to  (e_0,*[^e_0])*(e_0),  which  may  be  proved  by  induction  on  the  length 
of  a  sequence.  Writing  each  sequence  in  this  new  form  is  profitable  because  e_0  is  usually 
much  more  probable  than  the  other  events.  Therefore  it’s  faster  to  use  a  geometric 
distribution  to  find  the  lengths  of  the  e_0  runs,  and  to  only  use  the  inverse  transformation 
method  to  determine  the  events  that  fall  in-between  these  runs. 

3.5  PERFORMANCE 

3.5.1  Create  Spectra 

This  function  generates  an  eight-band  IBSM  file  containing  atmospheric  and  sensor 
spectral  distributions.  The  bands  contain,  respectively,  (0)  atmospheric  path  transmission, 
(1)  atmospheric  path  radiance,  (2)  irradiance  at  the  target  for  a  direct  source,  (3)  diffuse 
downwelling  radiance  at  the  target,  (4)  target  reflectance,  (5)  background  reflectance,  (6) 
sensor  optical  chain  transmission,  and  (7)  detector  quantum  efficiency.  All  are  in  single 
row  format  of  "nwavl"  columns.  When  the  override  parameters  are  not  selected,  the  first 
four  distributions  (0-3)  are  computed  via  MODTRAN  in  reference  to  the  ASCII  path  data 
file,  the  following  two  (4-5)  from  the  ASCII  reflectance  file  data,  and  the  final  two  (6-7) 
from  the  ASCII  sensor  file  data.  In  all  cases,  the  spectral  distributions  are  computed 
through  a  linear  interpolation  of  the  data,  and  can  be  overridden  with  the  spectrally  flat 
override  parameters.  The  minimum  and  maximum  x  scaling  attributes  contain  the 
wavelength  minimum  and  maximum  in  meters.  The  minimum  y,  maximum  y,  center 
wavelengths,  and  spectral  bandwidths  attributes  should  be  ignored.  This  generated  IBSM 
file  is  of  a  form  suitable  for  input  to  the  Sensor  Petformance  function. 
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3.5.2  Create MTF 


This  function  generates  a  one-dimensional  Modulation  Transfer  Function  (MTF) 
based  on  a  real-valued  two-dimensional  spatial  spectrum  file  representing  the  MTF.  The 
generated  MTF  can  lie  along  either  the  X  or  Y  axes  of  the  input  data,  and  is  placed  as  a 
single  row  in  an  IBSM  output  file.  The  frequency  scale  (min/max)  can  be  either  the  same 
as  the  input  (automatic)  or  overridden  (fixed).  In  the  latter  case,  the  data  is  linearly 
interpolated  to  the  specified  grid.  This  generated  IBSM  file  is  of  a  form  suitable  for  input 
to  the  Sensor  Performance  function. 

3.5.3  Sensor  Performance 


This  function  computes  metrics  characterizing  sensor  performance  (including  SNR, 
GRD,  and  NIIRS)  based  on  a  set  of  sensor  parameters;  target,  background,  atmosphere, 
and  sensor  spectral  distributions  using  the  Create  Spectra  function;  and  x  and  y  axis  MTF 
functions  using  the  Create  MTF  function.  The  metric  is  written  to  an  ASCII  file  along  with 
an  arbitrary  independent  variable.  Depending  on  the  "update"  parameter,  the  file  is 
initialized  before  writing  or  the  data  is  appended  to  an  existing  file.  In  the  long  form  output 
data  format,  the  full  input  parameter  set  is  written  to  the  file  upon  initialization  (update  =  0), 
and  the  metric  data  is  augmented  with  additional  intermediate  computation  results.  The 
short  form  output  data  format  results  in  a  two  column  ASCII  file  containing  the 
independent  variable  and  computed  metric  for  each  update  execution,  which  is  useful  for 
generating  parametric  plot  data. 


The  Signal  to  Noise  Ratio  (SNR)  metric  refers  to  the  target  to  background  detected 
signal  difference  relative  to  the  noise  level.  It  is  computed  by: 


SNR  = 


LOW 


Nfj,  +  (jf  + 


12 


N„:. 


-|2 


(3-47) 


where 


=  max(iV,^„A^,,J 


(3-48) 
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i,.,(^)  =  [l  -P,„(^)]i'..(^.  J’,„)+P,„(A)r-^=ii^^  +  i,,j,.„(A)]  (3-52) 


V(^)  =  [>  -  Pi>,W]b„(k.T,^)  +  P«.(A)[  +  L^.,.W]  (3-53) 


and  Lbb  (^,T)  is  the  blackbody  function  given  in  Eq.  (3-35).  The  parameters  used  in  Eqs. 
(3-47)  to  (3-53)  are  defined  in  Table  3-5. 

The  SNR  computation  also  includes  saturation  effects  and  optionally  automatic  gain 
control.  Saturation  is  incorporated  by  limiting  Nhi  and  Nlow  to  within  the  minimum  to 
maximum  range.  Automatic  gain  control  is  performed  by  reducing  both  Nhi  and  Nlow  by 
a  factor  such  that  Nhi  equals  a  specified  AGC  level.  If  Nhi  is  less  than  this  AGC  level 
initially,  no  scaling  is  performed. 

The  Ground  Resolved  Distance  (GRD)  metric  is  computed  orthogonal  to  the  sensor 
line-of-sight  axis.  Initially,  it  is  computed  in  both  the  x  and  y  axes  by 


R  ‘2-Pr.R 
Gm  =max  — 

«/  / 


(3-54) 
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Table  3-5:  SNR  Performance  Parameter  Definitions 


Ldiffuse  (^) 

Edirect  (^) 

Ptgt  (^)>  Pbkg  (^) 
Ttgt.  Tbkg 
Ltgt  (^)>  Lbkg  (^) 
Lpath  (^) 

'^atm  (^) 

Xopt (^) 

£opt 

Topt 

Ti  (X) 

Lstray 

Td 

NtdI 

f/# 

Nx,  Ny 
Wx,  Wy 
Njiiin^  Nmax 

b 

On 

Jd 

Ntgb  Nbkg 


Diffuse  downwelling  radiance  (W/m^sr) 
Direct  source  irradiance  (W/m^) 

Source  zenith  angle  (deg) 

Target,  background  reflectance 
Target,  background  temperature  (Kelvin) 
Target,  background  radiance  (W/m^sr) 
Atmospheric  path  radiance  (W/m^sr) 
Atmospheric  path  transmission 
Optical  train  transmission 
Effective  optics  emissivity 
Effective  optics  temperature  (Kelvin) 
Detector  quantum  efficiency 
In-band  stray  radiance  (W/m^sr) 

Dwell  time  (msec) 

Time  -  delay  -  integrate  stages 
Optics  F-number 
Aggregation  in  x  and  y 
Detector  width  in  x  and  y  (pm) 

Quantizer  minimum  and  maximum 
Quantizer  digitization  levels 
Detector  rms  noise  (electrons) 

Detector  dark  current  density  (pA/cm^) 
Target,  background  signal  level  (electrons) 
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where  R  is  the  sensor  to  target  slant  range,  px  and  py  are  the  detector  center-to-center 
spacings  in  x  and  y,  f  is  the  optical  system  focal  length,  and  Ur  is  the  highest  spatial 
frequency  to  satisfy 


MTF^(u\>—  (3-55) 

The  final  GRD  is  given  as  the  geometric  mean  (square  root  of  the  product)  of  the  GRD  in  x 
and  y. 


Image  quality  is  specified  by  a  National  Image  Interpretability  Rating  Scale  (NIIRS) 
value.  This  is  computed  from  the  SNR  and  modulation  transfer  functions  in  x  and  y  by  the 
Generalized  Image  Quality  Equation  [29]; 

e  =  11.81-1-3.32  logio  (RER  /  GSD)  -  1.48H  -  (3-56) 

SNR 

where 


H=4Hjr, 

(3-57) 

RER  =  ^RER^RERy 

(3-58) 

GSD  =  ^GSDfiSD^ 

(3-59) 

mzxER,  =  1.0  to  3.0) 

H  = 

leR„(i-25) 

monotonic 

not  monotonic 

(3-60) 

RER,,  =  ER,_, {0.5)- ER,, (-0.5) 

(3-61) 

GSD„y-f 

(3-62) 

ER,, (^)  =  0.5+  ‘  f'"'” 

;rJo  w 

.  (27tM>N 

sin  - ’•  dw 

V  /  J 

(3-63) 
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The  parameters  used  in  Eqs.  (3-56)  to  (3-63)  are  given  in  Table  3-6.  NIIRS  is  generally 
specified  in  reflective  bands  for  a  target  reflectance  of  0.15  and  background  reflectance  of 
0.07,  and  in  thermal  bands  for  blackbody  radiance  with  a  target  temperature  of  302K  and 
background  temperature  of  300K. 

3.5.4  Exposure  Control 

This  function  performs  exposure  control  by  determining  the  pixel  dwell  time,  TDI 
stages,  and  aggregation  levels  which  place  the  maximum  of  the  target  and  background 
signal  levels  between  specified  minimum  and  maximum  levels.  For  the  signal 
computation,  sensor  spectral  distributions  using  the  Create  Spectra  function  are  used.  The 
signal  computation  is  given  by  Equations  (3-50)  to  (3-53). 

Initially,  the  signal  level  is  computed  for  a  user  defined  desired  dwell  time,  the 
maximum  allowable  TDI,  and  the  minimum  allowable  aggregation.  If  the  maximum  of  the 
target  and  background  signal  level  (called  initial  signal)  falls  between  the  minimum  to 
maximum  signal  window,  the  initial  sensor  configuration  values  are  unchanged. 

If  the  initial  signal  exceeds  the  maximum  (saturation  condition),  first  the  TDI  stages 
are  decremented  by  steps.  Once  the  adjusted  signal  level  falls  below  the  maximum,  the 
process  stops  with  a  new  number  of  TDI  stages  having  been  established.  If  the  saturation 
condition  cannot  be  resolved  with  the  minimum  allowable  TDI,  the  dwell  time  is  reduced  to 
place  the  adjusted  signal  to  the  maximum  level  with  the  minimum  allowable  TDI  stages. 

If  the  initial  signal  is  less  than  the  minimum  level  (low  signal  condition),  first  the 
aggregation  levels  are  incremented  (alternately  x  and  y)  by  steps.  Once  the  adjusted  signal 
level  falls  above  the  minimum,  the  process  stops  with  the  new  aggregation  levels  having 
been  established.  If  the  low  signal  condition  cannot  be  resolved  with  the  maximum 
aggregation,  the  dwell  time  is  increased  to  place  the  adjusted  signal  to  the  minimum  level 
with  the  maximum  aggregation. 

The  output  consists  of  the  sensor  configuration  values,  which  are  placed  in  a  five- 
band  floating-point  VIFF  file  (one  value  per  band)  in  the  following  order:  (0)  pixel  dwell 
time  in  msec,  (1)  TDI  stages,  (2)  X  aggregation,  (3)  Y  aggregation,  (4)  adjusted  signal 
level.  The  values  can  be  extracted  to  global  variables  using  the  Print  Value  function. 
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Table  3-6:  NIIRS  Performance  Parameter  Definitions 


Pixel  position 

Px,y 

Detector  pixel  pitch  (pm) 

Nx,y 

Aggregation  in  x  and  y 

w 

Spatial  frequency  (rad’^)  in  x  or  y 

f 

Optical  system  focal  length  (cm) 

D 

Aperture  diameter  (cm) 

^min 

Minimum  wavelength  (pm) 

ERx,y(D 

Line  response  function  in  x  and  y 

R 

Target  to  sensor  slant  range  (km) 

GSDx,y 

Ground  sample  distance  in  x  and  y 

RERx,y 

Relative  edge  response  in  x  and  y 

Hx,y 

Edge  height  overshoot  in  x  and  y 

SNR 

Signal  to  noise  ratio 

G 

SNR  gain  due  to  image  processing 

RER 

Relative  edge  response  (mean) 

H 

Edge  height  overshoot  (mean) 

GSD 

Ground  sample  distance  (mean) 

Q 

Image  quality  (NIIRS) 
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3.5.5  Image  Quality 


This  function  computes  image  quality  metrics  including  line  response  functions, 
GSD,  SNR,  and  NIIRS  from  an  image  containing  test  patterns  and/or  edge  features.  The 
user  specifies  the  locations,  in  angular  (mrad)  units,  of  horizontal  and  vertical  edges  as  well 
as  ideally  uniform  image  regions  for  noise  estimation. 

First,  the  function  extracts,  centers,  and  normalizes  the  vertical  and  horizontal  edge 
response  functions.  The  image  is  assumed  to  exhibit  the  proper  sampling  characteristics 
such  that  these  will  correspond  to  ERx,y  (Q  in  Eq.  (3-63)  where  ^  corresponds  to  image 
pixels.  Next,  the  mean  and  variance  over  designated  target  and  background  regions  is 
computed.  Finally,  the  image  quality  metrics  are  computed  via  Equations  (3-56)  through 
(3-62)  along  with 

SNR  =  -^  (3-64) 

^max 

where  Ap  is  the  difference  between  the  mean  signal  levels  and  Omax  is  the  standard 
deviation  associated  with  the  region  exhibiting  a  higher  mean  signal. 

Based  on  user  options,  the  function  output  consists  of  the  two  edge  response 
functions  in  two  column  ASCII  format  and  an  ASCII  metric  summary  file.  Spatial  scaling 
(rad  units)  attributes  must  exist  for  the  input  image,  and  the  specified  edge  and  region 
locations  must  fall  within  the  angular  range  of  the  image. 

3.6  PROCESSING 

3.6.1  Image  Interlace 

This  function  interlaces  up  to  four  images  in  a  user-specified  direction.  The 
direction  can  be  either  along  the  x  axis  (width),  y  axis  (height),  or  in  bands  (element).  The 
attribute  information  is  also  appropriately  adjusted,  although  care  must  be  taken  such  that 
the  input  attribute  information  is  appropriate.  For  example,  nonsensible  output  scaling 
attributes  will  occur  if  the  input  scaling  data  was  not  representative  of  interleaved  images. 
The  interlaced  images  must  be  of  the  same  size  and  data  type. 
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3.6.2  Aggregation 


This  function  sums  image  elements  in  width,  height,  and/or  bands  to  a  user 
specified  aggregation  level.  By  virtue  of  the  aggregation,  the  output  image  is  smaller  in 
each  dimension  by  the  aggregation  factor.  The  residual  image,  when  the  size  is  not  an 
integer  multiple  of  the  aggregation  factor,  is  tmncated.  If  the  projection  option  is  selected 
for  a  particular  dimension,  the  aggregation  factor  is  overridden  by  the  dimension  size. 

Note  that  overflows  may  occur  in  the  summation  process. 

3.7  UTILITIES 

3.7.1  VIFFtoIBSM 

This  function  adds  attribute  information  to  a  VIFF  file  to  create  an  DBSM  format 
output.  The  spatial  scaling  and  spectral  band  data  can  either  be  entered  from  the  screen  (up 
to  four  bands),  or  input  from  an  ASCII  file.  The  ASCII  file  format,  displayed  in  Table  3- 
7,  corresponds  to  that  generated  using  the  List  IBSM  File  Data  function. 

3.7.2  IBSM  Info 

This  function  generates  an  ASCII  file  containing  information  extracted  from  the 
IBSM  header  and  attributes.  This  includes  data  type,  image  size,  number  of  bands,  spatial 
scaling,  and  the  center  wavelength  and  spectral  bandwidth  for  each  band.  The  file  format, 
displayed  in  Table  3-7,  is  appropriate  for  input  to  the  VIFF  to  IBSM  utility. 

3.7.3  Transform 

This  function  performs  conversions  between  the  image  and  spatial  spectrum 
domains  through  a  discrete  Fourier  transform.  The  forward  and  reverse  transformations 
are  given  by  [30] 


G(m,v)  =  Y  t 

y=0  ife=0  ^  ^ 


+  m) 


I(x,y)  = -  Z  Z  G{u.,vAe 

N  N  ^'=0  ^  ' 

'x*  ' y 


(3-65) 

(3-66) 
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Table  3-7: 

IBSM  Header  File  Format 

File  Name: 

fname 

File  Type: 

IBSM  Khoros  image  file  format 

VIFF  Information: 

Data  Type: 

data-type 

X_Pixels: 

X _pixels 

Y_Pixels: 

y _pixels 

Bands: 

bands 

IBSM  Information: 

X_Min: 

xjmin 

X_Max: 

xjmax 

Y_Min: 

yjnin 

Y_Max: 

yjnax 

Wavelengths 

Bandwidths 

wavelength  (1) 

• 

bandwidth  (1) 

• 

wavelength  (bands)  bandwidth  (bands) 
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where  I(x,y)  is  an  image,  G(u,v)  is  the  corresponding  spatial  spectrum,  (x,y)  are  angular 
coordinates  (rad),  (y,v)  are  angular  spatial  frequency  coordinates  (rad'l),  and  Nx  and  Ny 
are  the  transform  sizes  in  x  and  y. 


For  computational  efficiency,  the  transforms  are  performed  using  a  Fast  Fourier 
Transform  (FFT)  algorithm.  The  use  of  this  algorithm  places  some  constraints  on  the 
transform  size  sampling.  The  transform  sizes  must  be  a  power  of  two  in  x  and  y.  If  the 
image  or  spectrum  to  be  transformed  is  not  an  exact  power  of  two  in  size,  a  larger 
transform  size  must  be  used.  The  input  image  or  spectrum  is  automatically  zero  padded  to 
this  size.  The  sampling  of  the  transform  is  also  automatically  computed  based  on  the 
sample  spacing  of  the  input.  Specifically, 


^min 


U  = 

max 


V  = 

min 


V  = 

'^max 


1 

2 

•^max  ^min 

1 

2 

*^max  •^min 

-N^ 

1 

2 

^max 

iVmin 

f^y 

-1 

J 

1 

2 

y max  y min 

(3-67) 

(3-68) 

(3-69) 

(3-70) 


for  the  forward  transform  case  where  the  image  scale  factors  (xmji},  Xi][jqx>  yniin>  ymsix) 
correspond  to  after  the  zero  padding  operation.  Note  that  the  transform  is  always 
essentially  centered,  which  is  accomplished  by  quadrant  swapping  prior  to  and  after  the 
FFT. 


Spatial  scaling  attributes  must  be  available  for  the  input  image,  although  the  image 
data  units  are  arbitrary.  Both  complex  and  real  input  and  output  data  types  are  supported. 
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3.7.4  Create  Plot  Data 


This  function  extracts  a  row  or  column  of  an  IBSM  file  and  outputs  it  along  with 
spatial  scaling  data  (extracted  from  the  scaling  attributes)  into  two  column  ASCII  format  for 
plotting. 

3.7.5  Create  Radiance  Image 

This  function  converts  a  ground  reflectance  and,  optionally,  temperature 
distribution  into  a  radiance  image,  including  both  reflected  and  radiated  components.  The 
reflected  component  uses  only  the  reflectance  image  input  (assumed  Lambertian)  with  direct 
and  diffuse  illumination  computed  via  MODTRAN  in  reference  to  the  ASCII  path  data  fde. 
This  file  can  be  generated  using  the  MODTRAN  Input  tool.  It  is  important  that  the 
wavelength  limits  specified  in  the  MODTRAN  path  data  file  cover  each  spectral  band 
specified  in  the  input  image  header.  The  direct,  diffuse,  or  both  illumination  values  from 
MODTRAN  can  be  overridden.  The  radiated  component  is  computed  using  both  the 
reflectance  and  temperature  (Kelvin  units)  images,  and  is  ignored  when  the  latter  does  not 
exist. 


Mathematically,  the  output  image  is  given  by 


Aq+AA/2  p 


I{x.y)  =  pix.y)  f 

J  K 


An-A2/2 


dX 


An+AA</2 


(3-71) 


■[l-p{x,y)]  jLBBl^J{x,y)\dX 


Xq-AX/2 


where  p(x,y)  is  the  reflectance  image,  T(x,y)  is  the  temperature  image,  Edkect  (X)  is  the 
directional  spectral  illumination,  (|)d  is  the  zenith  angle  of  the  directional  source,  Ldiffuse  (^) 
is  the  diffuse  downwelling  spectral  radiance,  and  Nk  is  the  spectral  bandwidth  for  the  band 

of  interest. 


Spectral  band  attributes  must  be  available  in  the  reflectance  input  image,  and  the 
data  must  be  in  reflectance  units.  The  temperature  input  file  data  must  be  in  units  of 
absolute  temperature  (Kelvin).  The  output  image  data  is  in  radiance  (W/m^sr)  units. 
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3.7.6  Ground/Angle  Transform 


This  function  performs  a  coordinate  transformation  of  a  radiance  image  between  the 
horizontal  ground  plane  (meter  units)  and  line-of-sight  angular  (radian  units)  coordinates. 
Only  the  image  scaling  attribute  values  are  altered  based  on  the  viewing  geometry.  The 
defined  azimuth  axis  is  the  non-squinted  axis  (normal  to  the  line-of-sight  and  vertical  axes). 
The  scaling  of  the  squinted  axis  is  based  on  a  linear  projection,  which  occurs  under  the 
assumption  that  the  ground  field-of-view  in  the  squinted  direction  is  much  smaller  than  the 
slant  range.  No  image  distortion  or  perspective  effects  are  modeled. 

If  the  input  image  is  scaled  in  meter  units  in  the  viewing  plane  (orthogonal  to  line- 
of-sight)  rather  than  the  horizontal  plane,  this  utility  can  be  used  with  a  zenith  angle  of  zero 
to  convert  to  angular  units. 

Spatial  scaling  attributes  in  meter  units  must  exist  for  the  input  image.  The  spectral 
band  attributes  and  image  data  are  arbitrary. 

3.7.7  Load  Variables 

This  function  loads  a  set  of  global  variable  definitions  from  an  ASCII  file.  The  file 
format  parallels  the  "Global  Variables"  input  form,  and  include  simple  definitions  of 
variables  as  constants,  expressions,  and  cormnents.  Comments  are  denoted  by  a  "#"  at  the 
start  of  a  line. 

3.7.8  XgraphPlot 

This  function  provides  a  user  interface  to  the  Xgraph  plotting  application,  and 
allows  up  to  four  plots  per  graph  with  preformatted  titles,  labels,  and  axes  ranges.  The 
data  is  input  as  separate,  two-column  ASCII  files.  Xgraph  is  executed  from  a  system  call. 
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4.0  USAGE  OVERVIEW 


In  this  section,  information  is  provided  in  an  effort  to  assist  the  user  in  developing 
and  executing  sensor  models  based  on  the  IBSM  toolbox.  Detailed  information  on  the  use 
of  Khoros  and,  specifically.  Cantata  is  available  elsewhere  [31]  and  is  not  repeated  here. 
However,  brief  descriptions  of  some  commonly  used  features  are  given. 

4.1  FILE  FORMATS 

Image  data,  as  well  as  other  data  at  times,  is  passed  between  IBSM  modules  and 
other  Cantata  glyphs  (Khoros  term  for  functional  module)  by  files  in  the  Khoros 
Visualization/Image  File  Format  (VlFF).  The  Khoros  Image  File  Format  is  organized  as 
1024  bytes  of  header  followed  by  map(s),  location  data  and  then  image  data.  The  header 
also  contains  information  identifying  the  header  format  used.  This  file  format  or  data 
structure  has  evolved  from  originally  supporting  only  images  to  supporting  multi¬ 
dimensional  data.  The  multidimensional  aspect  is  used  to  accommodate  multiband 
imagery. 

For  use  with  IBSM,  specific  header  information  relating  to  spatial  scaling  and 
spectral  band  characteristics  is  added  as  attributes  to  the  VlFF  files.  Attribute  data  is  merely 
arbitrary  appended  data  to  a  VlFF  file.  Many  of  the  IBSM  modules  require  this 
information  to  be  present  for  correct  operation.  This  extended  VlFF  file  format  is  referred 
to  as  IBSM  file  format.  The  VlFF  to  IBSM  Utility  is  used  to  add  this  attribute  information 
to  a  VlFF  file.  Since  the  attribute  data  is  essentially  extraneous  in  the  VlFF  file  format, 
IBSM  file  formats  are  compatible  with  all  the  standard  Khoros  functions  and,  in  most 
cases,  the  attribute  data  is  passed  through  undisturbed. 

There  are  four  basic  types  of  IBSM  attribute  information.  The  first  is  a  image  type 
flag:  IBSM_FLAG.  This  denotes  that  the  VlFF  file  is  of  IBSM  format  (IBSM_FLAG  = 
85).  The  second  is  spatial  scaling  information:  IBSM_SCALE.  This  contains  four  values 
which,  in  order,  specify  Xmin,  Xn,ax»  Ymin,  and  Ymax-  Xmin  specifies  the  first  column  of 
the  image,  and  Ymm  the  first  row.  The  third  is  spectral  band  information: 
IBSM_WAVELENGTH  and  IBSM_BANDWIDTH.  Each  is  an  array  of  values  containing 
the  spectral  band  centers  and  bandwidths,  respectively,  for  each  band  of  image  data  in 
meter  units.  The  fourth  is  data  compression  attributes:  IBSM_COMPRESSION_TYPE, 
IBSM_COMPRESSION_BAND_SIZE,  and  IBSM_ORIGINAL_SIZE.  Table  4-1 
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Table  4-1 :  IBSM_COMPRESSION_TYPE  Attribute  Values 


Value  Compression  Type 

1  JPEG  (8  bit) 

2  JPEG  (12  bit) 


3  VQ 

4  DPCM 


5 


Wavelet 
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summarizes  the  IBSM_COMPRESSION_TYPE  attribute  values. 
IBSM_COMPRESSION_BAND_SIZE  is  given  in  bytes  and  roSM_ORIGINAL_SIZE  in 
width  and  height  elements. 

Image  data  is  stored  as  part  of  the  VIFF  5-dimensional  data  object.  The  multiband 
IBSM  imagery  utilize  the  WIDTH  dimension  for  the  x-axis,  HEIGHT  dimension  for  the  y- 
axis,  and  ELEMENT  dimension  for  bands. 

Table  4-2  summarizes  the  assumed  input  and  output  attribute  data  (not  including 
compression)  for  correct  operation  of  each  IBSM  function.  Table  4-3  summarizes  the  input 
and  output  file  types  and  assumed  units  (if  applicable)  for  each. 

4 . 2  MODEL  GENERATION 

Cantata  is  a  graphical  programming  environment  for  generating,  in  the  case  of 
IBSM,  sensor  evaluation  models.  Each  of  the  IBSM  functions,  along  with  over  a  hundred 
standard  data  processing  and  visualization  programs  in  the  Khoros  system,  are  represented 
by  icons  (called  glyphs).  To  create  a  model,  the  user  selects  the  desired  glyphs  (and 
control  structures,  as  needed),  places  them  in  the  Cantata  workspace,  and  connects  them  to 
indicate  the  flow  of  the  model.  Such  a  model  is  called  a  workspace.  Workspaces  can  then 
be  executed,  saved,  and  restored  to  be  used  again  or  modified  later.  The  save  and  restore 
features  are  implemented  by  selecting  the  Workspace  button. 

4.3  DATA  AND  IMAGE  VIEWING 

Three  primary  data  and  image  viewing  capabilities  exist  under  Cantata,  all  of  which 
are  selected  under  the  Output  button.  Information  contained  in  ASCII  files  can  be 
displayed  using  the  ASCII  File  Viewer  cohtained  under  the  Information  heading.  Imagery 
can  be  displayed  in  a  variety  of  ways,  all  contained  under  the  Display  Image  heading. 
Finally,  data  contained  in  ASCII  or  VIFF  files  can  be  plotted  using  the  Xprism  2  Utility 
under  the  Plot  Data  heading  or  using  the  Xgraph  Plot  IBSM  utility.  Plots  are  placed  on  the 
display  in  a  default  mode,  but  can  be  altered  by  the  user  and  output  to  a  printer  or  file. 
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Table  4-2:  Assumed  IBSM  Attribute  Information  (Continued) 


Table  4-3:  Input  and  Output  File  Types  and  Units 
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Table  4-3:  Input  and  Output  File  Types  and  Units  (Continued) 
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4.4  STANDARD  KHOROS  TOOLBOXES 


In  addition  to  the  IBSM  functions  incorporated  under  Khoros,  the  user  can  utilize  a 
broad  range  of  standard  utilities  for  signal  and  image  processing  by  selecting  glyphs  from 
the  standard  toolboxes.  A  brief  overview  is  given  here. 

4.4.1  Arithmetic 

This  toolbox  contains  routines  for  a  range  of  mathematical  image  operations, 
including  unary  operations  (scaling,  offset,  absolute  value,  square  root,  etc.),  binary 
operations  (addition,  subtraction,  etc.),  logical  operations  (And,  Or,  etc.),  and  linear 
transforms. 

4.4.2  Control 

This  toolbox  contains  routines  for  looping  and  conditional  constructs.  The  control 
functions  are  useful  particularly  for  looping  (e.g.,  computing  sensor  performance  as  a 
function  of  some  parameter)  and  conditional  branches.  The  generic  command  allows  the 
user  to  insert  an  arbitrary  UNIX  command  line  program  anywhere  within  a  workspace. 
Loops  and  procedures  (groups  of  glyphs  combined  into  one)  contain  all  included  glyphs  as 
a  sub- workspace. 

4.4.3  Data  Manipulation 

This  toolbox  contains  a  range  of  image  processing  functions  including  analysis 
(statistics),  data  conversion,  transforms,  frequency  filters,  histograms,  geometric 
operations  (expanding,  shrinking,  transposing,  etc.),  and  subregion  operations  (cropping, 
inserting,  etc.) 

4.4.4  Geometry 

This  toolbox  contains  routines  for  operating  with  geometrical  constructs. 
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4.4.5  Image  Processing 


This  toolbox  contains  routines  for  data  compression,  geometrical  operators,  and 
image  filtering. 

4.4.6  Input/Output 

This  toolbox  contains  utilities  for  inserting  supplied  or  user  defined  image  and 
signal  data  files  into  a  workspace,  or  for  creating  images  or  signals,  and  converting  file 
formats.. 

4.4.7  Matrix 

This  toolbox  contains  routines  for  matrix  manipulation. 

4.4.8  Program  Utilities 

This  toolbox  contains  generic  interfaces,  and  command  and  comment  icons. 

4.4.9  Visualization 

This  toolbox  contains  routines  for  image  and  data  capture  and  display. 

4.5  GLOBAL  VARIABLES 

The  Global  Variables  button  is  used  to  define  variables  and  evaluate  expressions 
within  the  Cantata  environment.  Once  defined,  variables  can  be  used  in  the  place  of 
integer,  float,  and  double  arguments  (GUI  entries)  of  glyphs.  Also,  expressions  can  be 
defined  to  establish  interrelationships  between  variables  (and,  therefore,  model  parameters) 
which  can  be  dynamic  in  the  course  of  a  workspace  or  model  execution. 

For  example,  consider  a  sensor  performance  model  which  loops  over  some 
independent  variable.  It  is  conceivable  that  the  IBSM  modules  in  each  iteration  of  the  loop 
change  dependent  on  the  loop  variable.  This  relationship  is  established  by  representing  the 
IBSM  module  input  as  a  variable  and  defining  the  variable  by  some  functional  dependence 
to  the  loop  variable.  Valid  expressions  for  this  functional  dependence  include  variables. 
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standard  arithmetic  operators  and  logicals,  as  well  as  predefined  constants  and  functions. 
Any  string  of  alphanumeric  characters,  beginning  with  a  letter,  may  be  used  as  a  variable. 

In  addition  to  defining  functional  dependencies  between  global  variables  and 
defining  the  glyph  operation  in  terms  of  such  variables,  it  is  also  possible  to  define  a 
dependence  of  a  global  variable  based  on  the  result  of  a  glyph  operation.  This  can  be 
accomplished  using  the  Print  Value  utility  under  the  Data  Manip  /  Analysis  and  Information 
heading,  which  allows  the  user  to  assign  any  element  of  a  VIFF  file  to  a  global  variable. 

4.6  ON-LINE  DOCUMENTATION 

On-line  documentation  (or  help)  is  provided  by  clicking  on  help  buttons.  The  on¬ 
line  help  pages  that  are  displayed  by  a  particular  help  button  are  appropriate  to  the  location 
of  the  graphical  user  interface  on  which  the  help  button  is  found.  Since  the  Cantata  master 
form  is  at  the  highest  level  of  the  Cantata  GUI  hierarchy,  the  help  button  at  the  upper  right 
hand  comer  of  the  Cantata  master  form  displays  by  default  the  introduction  to  Cantata. 
Notice  the  button  at  the  upper  left  hand  comer  of  the  on-line  help  display  labeled,  "More 
Help  Pages".  This  submenu  button  allows  you  to  display  related  help  pages;  in  this  case, 
all  the  other  help  pages  for  Cantata,  organized  not  alphabetically  but  in  the  order  in  which 
they  appear  in  the  manual. 

Help  buttons  found  on  each  subform  will  display  help  pages  which  describe  the 
purpose  and  operation  of  that  particular  subform.  Again,  any  other  related  help  pages  may 
be  accessed  by  using  the  submenu  button  labeled,  "More  Help  Pages." 

The  on-line  documentation  actually  places  the  man  pages  for  a  particular  function 
into  a  user  screen.  In  addition  to  a  description  of  the  function  and  usage,  a  description  of 
the  model  inputs  is  given.  Since  each  function  is  executable  from  a  UNIX  command  line, 
each  input  is  associated  with  a  unique  command  line  switch.  Therefore,  the  various  inputs 
in  the  documentation  are  indicated  by  these  switches.  The  association  with  GUI 
descriptions  is  not  explicitly  given,  but  is  probably  obvious  in  most  cases. 
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Figure  5-1:  Workspace  Corresponding  to  Test  Case 
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Figure  5-2:  Sensor  OTF  Chain  Procedure 
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Figure  5-3:  MODTRAN  Input  Form 
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Figure  5-4:  Comparison  of  Test  Case  MTFs  (X-axis) 
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Y  Axis  MTF  Comparison 
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Figure  5-5:  Comparison  of  Test  Case  MTFs  (Y-axis) 


Table  5-2:  GIQE  Test  Metric  Comparison 


GIQE 

Documentation 

IBSM 

Example 

GSD 

3.8  inches 

3.79  inches 

RER 

0.57 

0.582 

H 

0.94 

0.938 

G 

2.8 

2.8 

Signal 

24,000  electrons 

15,700  electrons 

Noise 

312  rms  electrons 

235  rms  electrons 

SNR 

77 

66.6 

NIIRS 

7.6 

7.75 
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5.2  MTF  ANALYSIS 


The  second  example  is  an  illustration  of  the  individual  MTF  components  for  the 
GIQE  test  case.  The  workspace  is  depicted  in  Figure  5-6.  In  this  case,  there  are  parallel 
paths  for  generating  the  MTF  functions  (in  the  X  direction)  due  to  the  atmospheric  path, 
optics,  line-of-sight  motion,  and  detector.  Figure  5-7  provides  the  resulting  MTF  plots. 

5.3  PARAMETRIC  PERFORMANCE  LOOP 

The  third  example  is  an  illustration  of  computing  sensor  performance  (SNR  and 
NIIRS)  as  a  function  of  an  independent  variable,  specifically  ground  range.  Essentially, 
the  model  from  Section  5.1  is  inserted  inside  a  loop  which  adjusts  ground  range  from  0  to 
100  km  in  5  km  increments.  The  workspace  is  given  in  Figure  5-8.  In  addition,  a 
conditional  branch  exists  which  sets  an  update  parameter  for  the  Sensor  performance 
modules  such  that  they  initialize  the  output  data  files  on  the  first  iteration  and  append  them 
on  subsequent  iterations.  The  resulting  SNR  and  NIIRS  plots  versus  ground  range  are 
given  as  Figures  5-9  and  5-10.  These  differ  slightly  from  the  test  case  because  MTF 
compensation  was  not  included  in  the  loop. 

5.4  IMAGE  BASED  ANALYSIS 

The  fourth  example  represents  the  primary  intended  use  of  the  IBSM  capabilities. 
Using  the  same  hypothetical  sensor  description,  the  workspace  depicted  in  Figure  5-1 1 
simulates  an  output  image  from  a  high  resolution,  high  SNR  input  scene,  and  extracts  edge 
response  functions  from  a  test  border  around  the  image  to  compute  NIIRS. 

The  input  scene  is  appropriately  scaled  spatially  and  in  (apparent)  reflectivity  units. 
The  model  proceeds  as  follows.  First,  an  exposure  control  sequence  is  executed  to  avoid 
deteetor  saturation.  Next,  the  input  scene  is  converted  into  a  line-of-sight  radiance  image 
and  propagated  through  the  atmospheric  path.  The  image  is  then  transformed  to  the  spatial 
spectrum  domain  and  passed  through  a  series  of  OTFs.  After  transforming  back  to  the 
image  domain,  the  radiance  image  is  converted  to  a  focal  plane  irradiance  image,  then  a 
detected  photoelectron  image,  sampled  at  the  detector  pitch,  perturbed  with  random  noise, 
and  quantized.  Finally,  the  edge  response  functions  are  extracted  and  sensor  performance 
is  computed. 
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Figure  5-6:  Workspace  Corresponding  to  MTF  Analysis 
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Figure  5-7:  MTF  Analysis  Results  (X-axis) 
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Figure  5-8:  Workspace  Corresponding  to  Parametric  Performance  Loop 


Figure  5-9:  Signal-to-Noise  Ratio  Versus  Ground  Range 
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Figure  5-10;  Image  Quality  Versus  Ground  Range 
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Figure  5-11:  Workspace  for  Image  Based  Analysis 


The  input  scene  is  shown  in  Figure  5-12,  and  can  be  compared  to  the  simulated 
output  image  shown  in  Figure  5-13.  The  edge  response  functions  are  illustrated  in  Figure 
5-14,  and  the  sensor  performance  computations  are  summarized  in  Table  5-3. 

Utilizing  an  appropriate  input  scene  perspective  and  adjusting  the  ground  range  to 
100  km,  an  analogous  simulation  with  a  greater  amount  of  image  degradation  is  executed. 
The  results  are  shown  in  Figures  5-15  to  5-17  and  Table  5-4. 

5.5  DATA  LINK  MODELING 

The  final  example  is  an  illustration  of  utilization  of  the  data  link  toolbox  which  has 
been  added  to  IBSM.  Figure  5-18  provides  the  workspace  for  modeling  the  data  link 
summarized  in  Table  5-5,  consisting  of  JPEG  compression/decompression,  a  finite  buffer, 
and  channel  errors.  The  input  scene  is  shown  as  Figure  5-19.  This  image  exhibits  a  fairly 
nonuniform  amount  of  compression  over  the  image  since  the  image  is  much  more  detailed 
near  the  center.  The  bits/block  file,  shown  as  Figure  5-20,  illustrates  this  fact.  For  the 
chosen  parameters  and  image,  the  buffer  overflows  near  the  end  of  the  image  (4096th 
block),  as  shown  in  the  buffer  occupancy  plot  in  Figure  5-21.  The  output  image  (Figure 
5-22)  illustrates  terminal  failure  of  the  decompressor  at  this  point  due  to  the  lost 
information,  as  well  as  temporary  failure  for  channel  bit  errors  (recovery  at  the  restart 
codes). 
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Figure  5-12:  Input  Image  for  20  km  Ground  Range  Simulation 
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Figure  5-13:  Output  Image  for  20  km  Ground  Range  Simulation 
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Table  5-3:  Sensor  Performance  at  20  km  Ground  Range 


Relative  Edge  Response 
RERx  =  0.394 
RERy  =  0.371 
RER  =  0.383 

Edge  Height  Overshoot 
Hx  =  0.901 
Hy  =  0.893 
H  =  0.897 

Ground  Sample  Distance 
GSDx  =  3.793  inches 
GSDy  =  3.793  inches 
GSD  =  3.793  inches 

Signal  to  Noise  Ratio 

Tgtmean  =  73159.12 
Tgt  stdev  =  290.65 
Bkg  mean  =  44882.627 
Bkg  stdev  =  242.39 
SNR  =  97.29 

NIIRS  Rating  =  7.25 


80 


Figure  5-15:  Input  Image  for  100  km  Ground  Range  Simulation 
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Figure  5-16:  Output  Image  for  100  km  Ground  Range  Simulation 


Table  5-4:  Sensor  Performance  at  100  km  Ground  Range 


Relative  Edge  Response 
RERx  =  0.427 
RERy  =  0.437 
RER  =  0.432 

Edge  Height  Overshoot 
Hx  =  0.977 
Hy  =  1.019 
H  =  0.998 

Ground  Sample  Distance 
GSDx  =  13.741  inches 
GSDy  =  13.741  inches 
GSD  =  13.741  inches 

Signal  to  Noise  Ratio 

Tgt  mean  =  96869.36 
Tgt  stdev  =  352.00 
Bkgmean  =  89230.78 
Bkg  stdev  =  391.26 
SNR  =  21.70 

NIIRS  Rating  =  5.38 
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Table  5-5:  Data  Link  Modeling  Parameters 


Image  Compression 
Type: 

Image  Quality  Factor: 
Blocks  per  Restart  Code: 


Buffer 


Size 

Input  Rate 
Output  Rate 
Delay  to  Empty 


Channel  Errors 


Flip  Errors 
Insert  Errors 
Delete  Errors 
Burst  Errors 
Burst  Length 


JPEG  (12  bit) 
3 
8 


64  KByte 

1  MWord/sec 

2  Mbits/sec 
6.4  msec 


10  errors/sec 

None 

None 

5  errors/sec 
1000  bits  (variable) 
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Figure  5-18:  Workspace  Corresponding  to  Data  Link  Effects  Simulation 
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Figure  5-20:  Bits/Block  Image  for  JPEG  Compressed  Image 
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6.0  SOFTWARE  OVERVIEW 


As  depicted  in  Figure  6-1,  the  IBSM  software  consists  of  primarily  three  functional 
layers  of  software:  a  Khoros  executive  layer  containing  file  I/O  and  graphical  user 
interface,  a  Khoros  library  function  layer  containing  functional  interface  between  executive 
layer  and  the  analytical  code,  and  the  sensor  modeling  library  containing  the  analytical 
modeling  code.  Numerical  Recipes  and  MODTRAN  are  used  for  analytical  computations 
and  atmospheric  modeling. 

6.1  KHOROS  EXECUTIVE  LAYER 

Table  6-1  sununarizes  the  routines  contained  in  the  Khoros  executive  layer.  Each 
routine  is  a  UNIX  command  line  executable  function,  but  also  contains  the  appropriate 
code  to  run  within  Cantata.  This  additional  code  includes  the  graphical  user  interface  and 
on-line  documentation.  The  appendix  provides  the  on-line  documentation  (Man  Pages)  for 
all  toolbox  functions. 

6.2  KHOROS  LIBRARY  FUNCTION  LAYER 

The  routines  contained  in  the  Khoros  library  function  layer  are  associated  with  the 
executive  layer  functions  in  Table  6-1 .  These  routines  are  C  functions  which  are  called  by 
Khoros  executive  functions.  The  majority  of  these  functions  contain  the  software  which 
loops  over  the  multiple  input  image  bands  and  calls  specific  analytical  function  to 
appropriately  process  each  image  band.  In  the  case  of  the  parametric  functions,  these 
functions  contain  module-specific  analytical  code  as  well. 

6.3  SENSOR  MODELING  LIBRARY 

Table  6-2  summarizes  the  routines  contained  in  the  sensor  modeling  library.  These 
routines  are  C  functions  which  form  the  analytical  modeling  basis  of  the  IBSM  software. 
These  routines  are  called  by  the  Khoros  library  functions  to  perform  the  majority  of  the 
physical  modeling  described  in  Section  3. 
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Figure  6-1:  IBSM  Software  Structure 


Table  6-1 :  Khoros  Executive  Layer  Routines 


Toolbox 

Function  Name 

Executive  Routine 

Library  Function  Routine 

Atmosphere 

Modtran  Input 

atmodtran 

Path  Trans/Rad 

atpatrra 

katmos_path 

Turbulence  OTF 

atturbotf 

kotf_turb 

Aero-Optic  OTF 

ataerotf 

kotf_aero 

Optics 

Diffraction  OTF 

opdiffotf 

kotf_diffr 

Defocus  OTF 

opdefotf 

kotf_blur 

Jitter  OTF 

opjittotf 

kotf Jitter 

Drift  OTF 

opdrifotf 

kotf_drift 

Wavefront  OTF 

opwaveotf 

kotf_wave 

User  OTF  (1-D) 

opuserotf 

kotf_user 

User  OTF  (2-D) 

opusrotfZ 

kotf_user2d 

Radiometry 

opradiom 

kradiometry 

Detector 

Detector  Size  OTF 

desizeotf 

kotf_detsize 

TDIOTF 

detdiotf 

kotf_tdi 

GTE  OTF 

decteotf 

kotf_cte 

Diffusion  OTF 

dediffotf 

kotf_detdiff 

Detection 

dedetect 

kdetection 

Noise 

denoise 

knoise 

Quantization 

dequantiz 

kquantize 

Sampling 

desamplin 

ksampling 

Nonuniformity 

denonunif 

knonunit 

Data  Link 

JPEG  Compress 

dljpegcomp 

jpeg_compress 

JPEG  Decompress 

dljpegdecomp 

jpeg_decomp 

TSVQ  Compress/Decompress 

dltsvq 

tsvq_(de)compress 

Wavelet  Comp/Decomp 

dlwavelet 

wave 

Finite  Buffer 

dlbuffer 

buffer 

Channel  Errors 

dlchannel 

channel 

Performance 

Create  Spectra 

pecrespec 

kspectra 

Create  MTF 

pecremtf 

kmtf 

Sensor  Performance 

pesensper 

kperf 

Exposure  Control 

peexpcd 

kexp_cntrl 

Image  Quality 

peimgqual 

kqimage 

Processing 

Image  Interlace 

printerlace 

Aggregation 

praggreg 

Utilities 

VIFFtoIBSM 

viffZibsm 

List  IBSM  File  Data 

ibsminfo 

Transform 

utxform 

ktransform 

Create  Plot  Data 

utcreplot 

Create  Radiance  Image 

utcreradi 

kradimg 

Ground/Angle  Transform 

utgrdangl 

kgndang 

Image  Interlace 

utinterlace 

Load  Variables 

utloadvars 

Xgraph  Plot 

utxgraph 
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Table  6-2:  Sensor  Modeling  Library  Routines 


Group 

Name 

Function 

Atmosphere 

atmos_path 

Apply  atmospheric  path  transmission  and  radiance  to  image 

otf_turb 

Generate  and  apply  path  turbulence  OTF  to  spatial  spectrum 

otf_aero 

Generate  and  apply  aero-optic  boundary  layer  OTF  to  spatial  spectrum 

atmos_rO 

Compute  atmospheric  path  correlation  diameter 

atmos_cndata 

Compute  index  structure  profile  along  path 

modtran_path 

Compute  atmospheric  transmission  and  path  radiance 

modtran_rad 

Compute  solar  or  lunar  direction  and  diffuse  radiance 

Optics 

otf_diffr 

Generate  and  apply  aperture  diffraction  OTF  to  spatial  spectrum 

otf_blur 

Generate  and  apply  gaussian  blur  OTF  to  spatial  spectrum 

otfjitter 

Generate  and  apply  line-of-sight  jitter  OTF  to  spatial  spectrum 

otf_drift 

Generate  and  apply  line-of-sight  drift  OTF  to  spatial  spectrum 

otf  wave 

Generate  and  apply  wavefront  error  OTF  to  spatial  spectrum 

otf_user 

Generate  and  apply  1-D  user  defined  OTF  to  spatial  spectrum 

otf_user2d 

Generate  and  apply  2-D  user  defined  OTF  to  spatial  spectrum 

radiometry 

Perform  optical  system  radiometric  operations 

Detector 

detection 

Perform  detector  photoelectronic  conversion 

gauss_noise 

Apply  gaussian  distributed  noise  to  image 

shot_noise 

Apply  signal  dependent,  Poisson  distributed  noise  to  image 

otf_detsize 

Generate  and  apply  detector  footprint  OTF  to  spatial  spectrum 

otf_detdiff 

Generate  and  apply  detector  diffusion  OTF  to  spatial  spectrum 

otf_tdi 

Generate  and  apply  time-delay-integrate  OTF  to  spatial  spectmm 

otf_cte 

Generate  and  apply  charge  transfer  inefficiency  to  spatial  spectrum 

sample 

Sample  image  based  on  detector  pitch 

quantize 

Quantize  image  based  on  scalar  truncation 

nonunif 

Apply  effects  of  gain  and/or  offset  nonuniformity 

Utilities 

transform 

Convert  between  image  and  spatial  spectrum  domains 

Miscellaneous 

copyhdr 

Copy  header  structure 

blackbody 

Compute  blackbody  spectral  radiance 

linint 

Perform  linear  interpolation 

lect 

Compute  rect  function 

sine 

Compute  sine  function 

tri 

Compute  tri  function 
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Appendix  -  Man  Pages 


Object  Information 

Category:  IBSM  Subcategory:  Atmosphere 

Operator:  Aero-Optic  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  aero-optic  OTF  to  a  spatial  spectrum  for  a  specified  turbulent 
boundary  layer.  Spatial  scaling  (rad-1  units)  and  band  attributes  must  be  available  for  the  input  image, 
and  the  data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced  from  an  image  using  the 
"Transform"  tool. 


Command  Line  Arguments 

Usage  for  ataerotf:  OTF  due  to  boundary  layer  aero-optic  turbulence 
%  ataerotf 

-i  (infile)  Name  of  input  IBSM  spatial  spectrum  file 

-o  (out file)  Name  of  output  IBSM  spatial  spectrum  file 

-vel  (float)  Sensor  platform  velocity  in  m/s 
vel  >=  0.0 

-rho  (float)  Air  density  in  kg/m3  at  sensor  altitude 
rho  >=  0.0 

-tbl  (float)  Boundary  layer  thic'kness  in  m  at  sensor  aperture 
tbl  >=0.0 


A.2.  atmodtran  —  Generate  a  MODTRAN  Path  Data  File 


Object  Information 

Category:  IBSM 

Operator:  MODTRAN  Input 

In  Cantata:  Yes 


Subcategory:  Atmosphere 

Object  Type:  kroutine 


Man  Pages 


IBSM  Manual  -  Chapter 


Description 

This  function  provides  a  user  interface  for  generating  a  MODTRAN  path  data  file  as  an  input  to  IBSM 
functions  that  execute  MODTRAN,  including  "Path  Trans/Rad",  "Create  Spectra",  and  "Create  Radi¬ 
ance  Image".  The  output  file  is  in  ASCII  format. 


Command  Line  Arguments 


Usage  for  atmodtran:  Generate  a  MODTRAN  Path  Data  File 
%  atmodtran 


-o 

-model 


(outfile)  Name  of  output  ASCII  MODTRAN  path  data  file 
(cycle)  Model  atmosphere 
0  "Tropical", 

1  "Midlatitude  Summer", 

2  "Midlatitude  Winter", 

3  "Subartic  Summer", 

4  "Subartic  Winter", 

5  "1976  US  Standard" 


-iday  (integer)  Day  of  year  in  days  starting  Jan.  1 

1  <=  iday  <=  365 


-imult 

(cycle) 

Multiple  scattering  switch 

0 

"Off" , 

1 

"On" 

-isourc 

(cycle) 

Celestial  source 

0 

"Sun" , 

1 

-Moon" 

“ihaze 


(cycle)  Aerosol  haze  model 
0  " None " , 

1  "Rural  23  km", 

2  "Rural  5  km", 

3  "Navy  Maritime", 

4  "Maritime  23  km" , 

5  "Urban  5  km" , 

6  "Troposheric  50  km", 

7  "Fogl  0.2  km", 

8  "Fog2  0 . 5  km" , 

9  "Desert" 


-anglem 

-icld 


(float)  Moon  phase  auigle  in  degrees 
unbounded 

(cycle)  Cloud  model 
0  "None", 

1  "Cumulus", 

2  "Altostratus" , 

3  "Stratus", 

4  "Stratocxomulus" , 

5  "Nimbostratus" , 

6  "Drizzle", 

7  "Light  Rain", 

8  "Moderate  Rain", 

9  "Heavy  Rain", 

10  "Extreme  Rain", 

11  "Standard  Cirrus", 

12  "Subvisual  Cirrus", 

13  "NOAA  Cirrus" 
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-parml  (float)  Source  azimuth  angle  in  degrees  (behind  sensor  =  180) 
unbounded 

-vis  (float)  Visibility  in  km  (0.0  =  default  for  haze  model) 

vis  >=  0.0 

-parni2  (float)  Source  zenith  angle  in  degrees  (vertical  =  0) 
unbounded 

-rainrt  (float)  Rain  rate  in  mm/hr 
rainrt  >=  0.0 

-wavll  (float)  Minimum  wavelength  in  microns 
unbounded 

-wavl2  (float)  Maximum  wavelength  in  microns 
unbounded 

-tbound  (float)  Surface  temperature  in  K 
tbound  >=  0.0 

-idv  (integer)  Spectral  sampling  in  cm-1 

1  <=  idv  <=  2.14748e+09 

-ifwhm  (integer)  Spectral  resolution  in  cm-1 

2  <=  ifwhm  <=  2.14748e+09 

-hi  (float)  Sensor  altitude  in  km  above  ground  level 

unbounded 

-gndalt  (float)  Ground  altitude  in  km  above  sea  level 
unbounded 

-h2  (float)  Target  altitude  in  km  above  ground  level 

unbounded 

-angle  (float)  Target  zenith  angle  in  degrees  (180  =  downlooking) 
unbounded 

[-i]  (infile)  trigger  input  to  sync  execution 

(default  =  {none}) 

[-salbf]  (float)  Surface  albedo 

(unbounded,  default  =  0.2) 

[-salbp]  (cycle)  Albedo  model  index 

0  “ Snow  Cover " , 

1  "Forest", 

2  “Farm", 

3  "Desert", 

4  "Ocean" , 

5  "Cloud  Deck" 

[default;  0  "Snow  Cover"] 


A.3.  atpatrra  —  Atmospheric  path  transmission  and  radiance 


Object  Information 

Category:  IBSM  Subcategory:  Atmosphere 

Operator:  Path  Trans/Rad  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  applies  the  effects  of  atmospheric  path  transmission  and  radiance  (scattered  and  radiated) 
globally  to  a  radiance  image.  The  function  executes  MODTRAN  in  conjunction  with  path  data  con¬ 
tained  in  the  MODTRAN  path  data  file.  This  file  can  be  generated  using  the  "MODTRAN  Input"  tool. 
It  is  important  that  the  wavelength  limits  specified  in  the  MODTRAN  path  data  file  cover  each  spectral 
band  specified  in  the  input  image  attributes.  The  spatial  scaling  attributes  in  the  input  image  are 
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ignored. 


Command  Line  Arguments 

Usage  for  atpatrra;  Atmospheric  path  transmission  and  radiance 
%  atpatrra 

-i  (infile)  Ncune  input  IBSM  radicuice  image  file 

-o  (out file)  Name  of  output  IBSM  radiance  image  file 

-ipath  (infile)  Name  of  input  ASCII  MODTRAN  path  data  file 

[-workdir]  (string)  Optional  working  directory  for  MODTRAN 
(default  =  {none}) 


A.4.  atturbotf  —  OTF  due  to  long  path  atmospheric  turbulence 


Object  Information 

Category:  IBSM  Subcategory:  Atmosphere 

Operator  Turbulence  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  a  long  path  turbulence  OTF  to  a  spatial  spectrum.  The  turbulence 
magnitude  is  specified  by  the  imaging  geometry,  vertical  index  structure  parameter  model,  and  refer¬ 
ence  index  stmcture  parameter  at  h=lm.  Spatial  scaling  (rad-1  units)  and  band  attributes  must  be  avail¬ 
able  for  the  input  image,  and  the  data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced 
from  an  image  using  the  "Transform"  tool. 


Command  Line  Ai^pments 

Usage  for  atturbotf:  OTF  due  to  long  path  atmospheric  turbulence 
%  atturbotf 

-i  (infile)  Name  of  input  IBSM  spatial  spectrum  file 

-o  (out file)  Name  of  the  output  IBSM  spatial  spectrxam  file 

-hi  (float)  Target  altitude  in  km 

unbounded 

-h2  (float)  Sensor  altitude  in  3an 

unbounded 

-zenith  (float)  Target  zenith  angle  (downlooking  =  180  deg) 

unbounded 

-model  (cycle)  Vertical  index  structure  parameter  model 

0  "Sirtple", 

1  "Hufnagel", 

2  "Empirical" 

-V  (float)  Upper  atmosphere  wind  speed  in  m/s  (Hufnagel  model) 

unbounded 

-cn2_lm  (float)  Reference  index  structure  parameter  at  h  =  Im 

cn2_lm  >=  0.0 

-vrel  (float)  Sensor  velocity  relative  to  the  atmosphere  in  m/s 
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unbounded 

-td  (float)  Integration  time  in  msec 

unbounded 

-d  (float)  Aperture  diameter  in  cm 

d  >=  0.0 

[-rOatm]  (float)  Path  correlation  diameter  in  cm  (override) 
(rOatm  >=  0.0,  default  =  0) 


B.  Data  Link  Man  Pages 


B.l.  dlchannel  —  Simulates  data  channel  errors 


Object  Information 

Category:  IBSM  Subcategory:  Data  Link 

Operator.  Channel  Errors  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  inputs  an  IBSM  format  image  and  simulates  passing  it  through  a  noisy  channel.  It  then 
outputs  the  resulting  corrupted  image. 


Command  Line  Arguments 


Usage  for  dlchannel:  Simulates  data  channel  errors 
%  dlchannel 

-i  (infile)  Input  file  that  will  be  corrupted  by  the  channel  simulation 

-o  (outfile)  Output  file  which  contains  the  corrupted  data  by  the  channel  simulatd 

-pf  (double)  The  probability  that  a  zero  will  be  changed  to  a  one,  or  a  one  will  1 

0  <=  pf  <=  1 

-pi  (double)  The  probability  that  a  bit  inserted  at  any  particular  position.  This 

0  <=  pi  <=  1 

..pd  (double)  The  probability  that  any  particular  bit  will  be  deleted. 

0  <=  pd  <=  1 

-pb  (double)  The  probability  that  a  burst  error  will  start  at  any  particular  bit. 

0  <=  pb  <=  1 

.nib  (double)  The  mean  length  of  the  error  bursts.  Together  with  the  Distribution  c 

mb  >=  0.0 

-db  (cycle)  Determines  whether  the  burst  errors  have  constant  or  variable  length 

0  "Constant  Length" , 

1  "Variable  Length" 


(cycle)  If  this  is  "zero"  or  "one"  then  all  of  the  bits  in  a  burst  error  wil] 
0  "Zeros", 

1  "Ones", 

2  "Zeros /Ones* 
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-hdrlen  (integer)  Defines  the  number  of  bits  at  the  beginning  of  each  band  that  are  gua 
hdrlen  >=  0 


B.2.  dlbuffer  —  Simulates  the  effects  of  buffering  on  input  data 


Object  Infonnation 

Category:  IBSM  Subcategory:  Data  Link 

Operator  Finite  Buffer  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  simulates  the  buffering  of  an  IBSM  format  image.  It  inputs  an  image  and  its  correspond¬ 
ing  bits/block  image  (see  the  dljpegcomp  man  page  for  more  information).  The  latter  is  produced  auto¬ 
matically  by  IBSM  compressor  glyphs,  or  it  may  be  produced  manually  for  non-compressed  images.  Its 
outputs  consist  of  the  buffered  data,  which  may  be  different  from  the  the  original  data  if  overflows  and 
underflows  occurred.  It  also  produces  a  graphable  data  file  that  represents  how  many  bits  were  in  the 
buffer  at  each  time  step. 


Command  Line  Arguments 

Usage  for  dlbuffer:  Simulates  the  effects  of  buffering  on  input  data 
%  dlbuffer 

-i  (infile)  Input  IBSM  compressed  file  to  buffer 

-ibits  (infile)  Input  IBSM  bits/block  file  from  conpression  output 

-p_in  (double)  The  time  between  element  arrivals. 

p_in  >  0.0 

•"P— out  (doxible)  The  time  between  element  departures, 

p_out  >  0.0 

-mode  (cycle)  Configures  the  buffer  to  operate  as  if  it  was  either  on  the  input  or 

0  "Input  Bit /Output  Block", 

1  "Input  Block/Output  Bit" 

-buf_s  (integer)  The  total  amount  of  bits  the  buffer  cam  hold. 

buf_s  >=  0 

-sim_un  (cycle)  Turns  on  or  off  the  simulation  of  buffering  underflow, 

0  "Off", 

1  "On* 

-sim_ov  (cycle)  Turns  on  or  off  the  simulation  of  buffering  overflow. 

0  "Off", 

1  "On" 

-delay  (double)  The  number  of  seconds  that  the  buffer  is  allowed  to  fill  before  the  i 
delay  >=  0.0 

-ilev  (integer)  In  effect  the  buffer  is  filled  with  this  many  dummy  bits  before  it  st 

ilev  >=  0 

-hdr_len  (integer)  Nxamber  of  bits  protected  from  buffer  corruption 
hdr„len  >=  0 
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[-0]  (outfile)  Resulting  output  data  file  after  buffering 

(default  =  {none}) 

[-obof]  (outfile)  Holds  a  one-dimensional  vector  filled  with  the  number  of  bits  in  the 
(default  =  (none)) 

[-step_sz]  (double)  If  secified,  the  buffer  occupancy  file  is  updated  once  every  step  siz 
(unbounded,  default  =  -1) 


B.3.  dljpegcomp  —  Performs  NITFS  JPEG  Compression 


Object  Information 

Category:  IBSM  Subcategory:  Data  Link 

Operator:  JPEG  Compress  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  inputs  an  IBSM  format  image,  an  image  quality  factor,  and  a  restart  interval  length.  It 
compresses  each  band  separately  and  returns  both  the  compressed  data  and,  optionally,  a  bits/block 
image. 

Glossary 

Bits/Block  Image  (BB  image):  This  is  an  IBSM  term.  The  BB  image  for  some  arbitrary  image  T  has 
one  element  for  each  MCU  in  T.  This  element  contains  the  number  of  bits  in  the  compressed  represen¬ 
tation  of  that  MCU. 

The  BB  image  is  derived  by  counting  the  number  of  bits  in  each  restart  interval  and  averaging  this  over 
the  number  of  MCU’s  in  a  restart  interval.  The  BB  image  is  set  equal  to  this  average  for  the  entire 
restart  interval.  Therefore  the  shorter  the  restart  interval  length,  the  more  accurate  the  BB  image.  And 
(nearly)  perfectly  accurate  BB  images  are  produced  if  and  only  if  the  restart  interval  is  one  MCU  long. 

Image  (Quality  Factor  (IQF):  An  NIC  and  IBSM  term  which  denotes  a  number  from  1  to  5  that  controls 
the  accuracy  with  which  the  image  is  reconstructed.  Lower  values  mean  that  fewer  bits  are  used  in  the 
compressed  file.  Unfortunately  they  also  mean  that  the  reconstructed  image  is  cruder.  This  range  of  val¬ 
ues  (i.e.  1  to  5)  was  set  by  p.  60  of  MIL-STD-188-198A  (and  also  page  40  of  MIL-STD-2500 A) ,  but  is 
not  part  of  the  JPEG  standard  itself. 

JPEG:  JPEG  stands  for  Joint  Photographic  Experts  Group,  the  name  of  the  committee  that  came  up 
with  a  standardized  method  for  image  compression.  More  information  about  it  can  be  found  at: 
http://www.cis.ohio-state.edu/hypertext/faq/usenet/jpeg-faq/top.html  from  the  book:  W.  Peenebaker,  J. 
Mitchell,  JPEG_Still_Image_Data_Compression_  Standard,  Van  Nostrand  Reinhold,  N.Y.,  1993  or 
from  MIL-STD-188-198A. 

MCU:  MCU  is  a  JPEG  term  which  stands  for  Minimum  Coded  Unit.  To  define  it,  start  with  any  image 
and  divide  it  into  8x8  blocks  which  cover  the  whole  image  but  do  not  overlap.  Then  each  of  these 
blocks  is  an  MCU. 

NJC:  NJC  stands  for  the  NTTFS  JPEG  Compression  Software  Library.  These  are  some  routines  that 
were  written  by  the  Optivision  Co.  (1477  Drew  Ave..  Suite  102.  Davis,  CA  95616)  to  provide  an 
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example  of  software  that  implements  MIL-STD-188-198.  Release  2.2i  is  used  here.  This  library  was 
chosen  over  its  competitors  partly  because  it  may  eventually  provide  unique  features  (such  as  the  use  of 
forward  error  correcting  codes  for  data  protection)  that  may  be  of  interest  to  users  of  IBSM. 

Restart  Interval:  This  is  a  JPEG  term  which  denotes  a  collection  of  one  or  more  consecutive  MCU’s. 

After  the  compressed  data  for  each  restart  interval,  the  compressor  transmits  a  code.  This  may  allow  the 
receiver  to  regain  synchronization  when  a  bit  error  occurred  in  the  compressed  data  for  that  restart  inter¬ 
val.  In  other  words,  chaimel  errors  should  be  confined  to  the  restart  interval  in  which  they  occurred.  So 
for  optimal  error  protection  we  want  the  restart  intervals  as  short  as  possible.  The  tradeoff  is  that  each 
restart  code  takes  about  20  bits  to  transmit,  so  the  shorter  the  restart  interval  the  more  bits  are  spent  on 
overhead.  To  comply  with  paragraph  5.2.3.3.1.4  of  MIL-STD-188-198A,  an  image’s  restart  interval 
must  be  no  larger  than  the  number  of  MCU’s  it  takes  to  span  the  width  of  that  image. 

Command  Line  Arguments 

Usage  for  dljpegcoir^):  Performs  NITFS  JPEG  Conpression 
%  dljpegcomp 

-i  (infile)  The  name  of  the  file  that  will  be  compressed 

-ocf  (out file)  Special  format  IBSM  file  that  contains  conpressed  JPEG  data 

-iqf  (integer)  Image  quality  factor  where  5  is  the  highest  quality  and  1  is  the  lowe 

1  <=  iqf  <=  5 

-rilen  (integer)  MCU=Minimum  Coded  Unit  (One  8x8  block):  The  number  of  8x8  image  block 

rilen  >  0 

[-obf]  (out file)  An  integer  image  each  of  whose  elements  contain  the  number  of  bits  pe 

(default  =  {none}) 


B.4.  dljpegdecomp  —  JPEG  decompression  of  the  input  image 


Object  Information 

Category:  IBSM  Subcategory:  Data  Link 

Operator  JPEG  Decompress  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  inputs  an  IBSM  format  image  created  with  the  dljpegcomp  function  with  possibly  post¬ 
processing  by,  for  example,  the  dlbuffer  or  dlchannel  functions.  It  then  creates  the  corresponding 
decompressed  image.  See  the  dljpegcomp  man  page  for  more  information. 


Command  Line  Arguments 

Usage  for  dljpegdecomp:  JPEG  decompression  of  the  input  image 
%  dljpegdecomp 

-i  (infile)  IBSM  formatted  file  that  was  compressed  with  the  JPEG  Compression  gly 

-o  (out file)  Decompressed  JPEG  image  file 
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B.5.  dltsvq  —  TSVQ  -  Tree-Structured  Vector  Quantization  Compression  and  Decompres¬ 
sion 


Object  Information 

Category:  IBSM  Subcategory:  Data  Link 

Operator:  TSVQ  Comp/Decomp  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  has  two  modes.  The  first  reads  an  input  image  and  a  TSVQ  codebook.  It  compresses  the 
input  using  the  codebook  and  produces  both  an  image  of  indices  and  a  bits/(unit  time)  image  (suitable 
for  the  buffer  glyph).  Also,  the  restart  interval  setting  determines  the  distance  between  consecutive 
restart  markers  in  the  image  of  indices,  unless  it’s  zero  in  which  case  restart  markers  are  not  used. 

The  second  mode  reads  an  image  of  indices  and  outputs  thedecompressed  image. 


Command  Line  Arguments 


Usage  for  dltsvq:  TSVQ  -  Tree- Structured  Vector  Quantization  Conqpression  and  Deconpression 
%  dltsvq 

-i  (infile)  File  name  of  the  original  image  or  corr5)ressed  image  depending  on  the 

-o  (outfile)  File  name  of  the  output  deconpressed  image  or  output  compressed  image 

-dtype  (integer  toggle)  Input  data  type  that  defines  whether  the  input /original  file  : 

1  (8  Bit) , 

2  (12  Bit), 
or  3  (Float) 


-mode  (integer  toggle)  Selects  the  mode  of  operation  as  either  compress  or  decompress 

1  (Compress), 
or  2  ( Decompress ) 


-ri_len  (integer)  After  every  ri^len  consecutive  indices  the  compressor  will  insert  a  i 
ri_len  >=  0 

[-bits]  (outfile)  Optional  output  file  that  represents  a  mesh  whose  interstices  are  8x^ 
(default  =  (none)) 

[-stem]  (string)  Codebook  file  names  stem  or  root  where  the  actual  file  names  are  ster 
(default  =  {none}) 


B.6.  dlwavelet  —  Wavlet  Compression  and  Decompression 


Object  Information 

Category:  IBSM  Subcategory:  Data  Link 

Operator:  Wavelet  Comp/Decomp  Object  Type:  script 

In  Cantata:  Yes 
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Description 

Wavelet  coding  tries  to  approximate  the  compression  performance  of  very  high-dimensional  vector 
quantization  at  a  reasonable  complexity.  The  wavelet  decoder  takes  the  indices  that  the  VQ  produced 
and  reconstructs  the  image  in  the  transformed  coordinate  system.  Finally,  it  returns  the  image  to  the 
normal  coordinate  system  by  multiplying  by  the  inverse  wavelet  transform. 


Command  Line  Arguments 


Usage  for  dlwavelet:  Wavlet  Conpression  and  Deconpression 
%  dlwavelet 

-i  (infile)  If  mode  is  "c"  this  should  contain  a  float-format  input  image  whose  d 

-o  (outfile)  If  nc  is  1,  this  is  the  name  of  a  file  that  will  contain  wavelet  tran 


[-stem] 

[-mode] 

[-ril] 

[-lev] 

[-nc] 

[default: 


(string)  The  prefix  of  two  codebook  file  names:  one  ends  in  "^shape . vif f  “ , 
(default  =  {none)) 

( string)  Wavelet  mode . 

(default  =  c) 

(integer)  The  restart  interval  length. 

(ril  >=  0,  default  =  0) 

(integer)  The  number  of  levels  in  the  wavelet  decomposition. 

(lev  >  0,  default  =  1) 

(integer  toggle)  Output  training  data 
1  (No) , 
or  2  (Yes) 

1] 


whi 


C.  Detector  Man  Pages 


C.l.  decteotf  —  OTF  due  to  charge  transfer  inefficiency 


Object  Information 

Category:  IBSM  Subcategory:  Detector 

Operator:  CTE  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  OTF  for  the  effects  of  charge  transfer  inefficiency  (CTE)  in 
charge-coupled  device  (CCD)  detector  readout.  The  per  transfer  CTE  is  assumed  identical  in  the  x  and 
y  directions.  Spatial  scaling  (rad-1  units)  and  band  attributes  must  be  available  for  the  input  image,  and 
the  data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced  from  an  image  using  the  "Trans¬ 
form"  tool. 
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Command  Line  Arguments 


Usage  for  decteotf:  OTF  due  to  charge  transfer  inefficiency 


%  decteotf 


-o 

”PX 


-py 

-f 

-ntx 

-nty 

-phase 

-eps 


(infile)  Name  of  input  IBSM  spatial  spectrum  file 

(outfile)  Name  of  output  IBSM  spatial  spectriim  file 

(float)  Detector  X  pitch  in  microns 
px  >=  0.0 

(float)  Detector  Y  pitch  in  microns 

py  0*0 

(float)  Sensor  focal  length  in  cm 
f  >=  0.0 

(integer)  Number  of  transfers  in  the  X  direction 
ntx  >  0 

(integer)  Number  of  transfers  in  the  Y  direction 
nty  >  0 

(integer)  N\imber  of  clock  phases 
phase  >  0 

(float)  Charge  transfer  efficiency  (per  phase) 
eps  >=  0.0 


C.2.  dedetect  —  Photoelectric  conversion  at  the  detector  plane 


Object  Information 

Category:  IBSM  Subcategory:  Detector 

Operator  Detection  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  converts  an  input  focal  plane  irradiance  image  into  a  noiseless  detected  photoelectron 
image  based  on  detector  parameters. 


Command  Line  Arguments 


Usage  for  dedetect:  Photoelectric  conversion  at  the  detector  plane 
%  dedetect 


-1 

-o 

-eta 

-ntdi 

-agx 

-agy 

-wx 

-wy 


(infile)  Name  of  input  IBSM  irradiance  image  file 
(outfile)  Name  of  output  IBSM  photoelectron  image  file 
(float)  Detector  mean  quantum  efficiency  over  band 
xanbounded 

(integer)  Number  of  TDI  stages  in  use 
ntdi  >  0 

(integer)  Aggregation  in  X  direction 
agx  >  0 

(integer)  Aggregation  in  Y  direction 
agy  >  0 

(float)  Detector  X  width  in  microns 
wx  >=  0.0 

(float)  Detector  Y  width  in  microns 
wy  >=  0.0 
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-td  (float) 

td  >=  0.0 
-jd  (float) 

jd  >=  0.0 
-sat  (float) 

sat  >  0.0 


Detector  dwell  (integration)  time  in  msec 
Dark  current  density  in  microamps/ cm2 

Saturation  photoelectron  level  (full  well)  in  electrons 


C.3.  dediffotf  —  OTF  due  to  CCD  charge  diffusion 


Object  Infonnation 

Category;  IBSM  Subcategory:  Detector 

Operator:  Diffusion  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  OTF  for  the  effects  of  carrier  diffusion  in  a  charge  coupled 
device  (CCD)  as  a  function  of  material  properties  and  dimensions.  Note  that  the  absorption  coefficient, 
at  least  in  the  case  of  silicon,  is  wavelength  dependent.  Spatial  scaling  (rad-1  units)  and  band  attributes 
must  be  available  for  the  input  image,  and  the  ^ta  must  correspond  to  a  spatial  spectrum,  which  can  be 
produced  from  an  image  using  the  "Transform"  tool. 


Command  Line  Arguments 

Usage  for  dediffotf:  OTF  due  to  CCD  charge  diffusion 
%  dediffotf 

-i  (infile)  Name  of  input  IBSM  spatial  spectrum  file 

-o  (outfile)  Name  of  the  output  IBSM  spatial  spectrum  file 

-alpha  (float)  CCD  xnaterial  absorption  coefficient  in  mm-1 
alpha  >=  0.0 

-aid  (float)  CCD  depletion  layer  width  in  mm 

aid  >=0.0 

-alo  (float)  Carrier  diffusion  length  in  mm 

alo  >=  0.0 

-f  (float)  Sensor  focal  length  in  cm 

f  >=  0.0 


C.4.  denoise  —  Addition  of  detector  and  electronic  noise 


Object  Information 

Subcategory:  Detector 

Category;  IBSM 

Operator:  Noise 

In  Cantata:  Yes 

Object  Type:  kroutine 
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Description 

This  function  adds  a  random  noise  process  to  a  photoelectron  image.  Two  types  of  noise  are  modeled. 
The  first  is  signal  dependent,  Poisson  distributed  shot  noise,  which  can  be  either  enabled  or  disabled. 
The  second  is  signal  independent,  additive  noise  that  is  assumed  zero  mean,  gaussian  distributed  with  a 
specified  root  mean  square  (RMS)  level.  The  noise  sources  are  assumed  both  jointly  and  pixel-to-  pixel 
independent. 


Command  Line  Arguments 

Usage  for  denoise:  Addition  of  detector  and  electronic  noise 
%  denoise 

_i  (infile)  Name  of  input  IBSM  photoelectron  image  file 

-o  (oUtfile)  Name  of  output  IBSM  photoelectron  image  file 

-shot  (cycle)  Shot  noise  selection  switch 

0  “No*, 

1  "Yes‘ 

-sigma  (float)  RMS  additive  noise  in  electrons 

sigma  >=  0.0 


C.5.  denonunif  —  Adds  effects  of  random  detector  nonuniformity  to  a  detected  image 


Object  Information 

Category:  IBSM  Subcategory:  Detector 

Operator  Nonuniformity  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  applies  the  effects  of  detector  nonuniformity  to  an  input  image.  The  level  of  nonunifor¬ 
mity  is  specified  separately  for  gain  (multiplicative)  and  offset  (additive)  portions.  Both  are  modeled  as 
gaussian  random  processes  with  the  latter  scaled  relative  to  a  user  specified  maximum  signal  level.  The 
user  can  also  specify  an  input  autocorrelation  file  which  characterizes  the  spatial  dependencies.  If  not 
specified,  the  nonuniformity  is  assumed  pixel-to-pixel  independent. 


Command  Line  Arguments 

Usage  for  denonunif:  Adds  effects  of  random  detector  nonuniformity  to  a  detected  image 


%  denonunif 


-i 

(infile) 

Input  IBSM  image  file  name 

-o 

(outfile) 

Output  IBSM  image  file  name 

-gain 

(float) 
gain  >= 

Percent  gain  uniformity 

0.0 

-offset 

(float) 
offset  > 

Percent  offset  nonuniformity 
=  0.0 

-nitiax 

(float) 

nmax  >= 

Meocimum  signal  level 

0.0 
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[-rho]  (infile)  Input  IBSM  autocorrelation  file  name 
(default  =  (none)) 


C.6.  dequantiz  —  Scalar  image  quantization 


Object  Infoimation 

Category;  IBSM  Subcategory:  Detector 

Operator:  Quantization  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  quantizes  a  detected  photoelectron  image  through  a  scalar  roundoff  process  to  2**b  digi¬ 
tization  levels  linearly  distributed  between  the  user-specified  minimum  and  maximum  levels. 


Command  Line  Arguments 

Usage  for  dequantiz:  Scalar  image  quantization 
%  dequantiz 

-i  (infile)  Name  of  input  IBSM  photoelectron  image  file 

-o  (outfile)  Name  of  output  IBSM  photoelectron  image  file 

-nmin  (float)  Minimiim  quantization  level  in  electrons 
nmin  >=  0.0 

-nmax  (float)  Maximum  quauitization  level  in  electrons 
nmax  >=  0.0 

-bits  (integer)  Digitization  levels  as  power  of  2 
bits  >  0 


C.7.  desamplin  —  Focal  plane  image  sampling 


Object  Information 

Category:  IBSM 

Operator:  Sampling 

In  Cantata:  Yes 


Subcategory:  Detector 

Object  Type:  kroutine 


Description 

This  function  performs  a  focal  plane  sampling  of  an  input  image  by  an  array  of  point  detectors,  assum¬ 
ing  detector  width  effects  have  been  previously  incorporated.  The  resampling  process  is  achieved  using 
a  bi-linear  interpolation  of  the  input  image  data  with  user  specified  absolute  alignment  of  the  sampling 
grid  to  the  input  image.  This  alignment  is  specified  by  way  of  the  fractional  (relative  to  sample  spac¬ 
ing)  offsets. 
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Command  Line  Arguments 

Usage  for  desamplin:  Focal  plane  image  sampling 
%  desamplin 

-i  (infile)  Name  of  input  IBSM  image  file 

-o  (outfile)  Name  of  output  IBSM  image  file 

-px  (float)  Detector  X  pitch  in  microns 

px  >=  0.0 

-py  (float)  Detector  Y  pitch  in  microns 

py  >=  0.0 

-agx  (integer)  Aggregation  in  X  direction 
agx  >  0 

-agy  (integer)  Aggregation  in  Y  direction 
agy  >  0 

-f  (float)  Sensor  focal  length  in  cm 

f  >=  0.0 

-xoff  (float)  Fractional  sampling  offset  in  X  direction 
unbounded 

-yoff  (float)  Fractional  sanpling  offset  in  Y  direction 
unbounded 


C.8.  desizeotf  —  OTF  due  to  finite  detector  width 


Object  Information 

Category:  IBSM  Subcategory:  Detector 

Operator:  Detector  Size  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  OTF  from  the  finite  detector  size  to  an  input  spatial  spectrum. 
The  detector  is  assumed  to  have  a  rectangular  form.  Spatial  scaling  (rad-1  units)  and  spectral  and  band 
attributes  must  be  available  for  the  input  image,  and  the  data  must  correspond  to  a  spatial  spectrum, 
which  can  be  produced  from  an  image  using  the  "Trartsform"  tool. 


Command  Line  Arguments 

Usage  for  desizeotf:  OTF  due  to  finite  detector  width 
%  desizeotf 

-i  (infile)  Name  of  input  IBSM  spatial  spectrum  file 

-o  (outfile)  Name  of  output  IBSM  spectrum  file 

-dx  (float)  Detector  X  width  in  microns 
dx  >=  0.0 

-dy  (float)  Detector  Y  width  in  microns 
dy  >=  0.0 

-px  (float)  X  detector  pitch  in  microns 
px  >  0.0 

-py  (float)  Y  detector  pitch  in  microns 
py  >  0.0 

-agx  (integer)  Aggregation  in  X  direction 


114 


Man  Pages 


IBSM  Manual  -  Chapter 


agx  >  0 

-agy  (integer)  Aggregation  in  Y  direction 
agy  >  0 

-f  (float)  Sensor  focal  length  in  cm 

f  >=  0.0 


C.9.  detdiotf  —  OFF  due  to  detector  time-delay-integration 


Object  Information 

Category:  IBSM  Subcategory;  Detector 

Operator:  TDl  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  OTF  due  to  the  detector  time-delay-integration  (TDI)  process  to 
an  input  spatial  spectrum.  The  TDI  can  occur  in  either  the  x  or  y  directions,  and  may  exhibit  a  mis¬ 
match  between  the  TDI  clocking  rate  and  the  Une-of-sight  (LOS)  motion  rate.  Even  for  a  perfect  rate 
match  (beta=1.0),  an  OTF  effect  occurs  due  to  the  LOS  drift  of  one  detector  phase.  Spatial  scaling 
(rad-1  units)  and  band  attributes  must  be  available  for  the  input  image,  and  the  data  must  correspond  to 
a  spatial  spectrum,  which  can  be  produced  from  an  image  using  the  "Transform"  tool. 


Command  Line  Arguments 

Usage  for  detdiotf:  OTF  due  to  detector  time-delay-integration 
%  detdiotf 

“i  (infile)  Name  of  input  IBSM  spatial  spectrum  file 

-o  (outfile)  Ncime  of  output  IBSM  spatial  spectrum  file 

-dir  (cycle)  Direction  of  TDI  (X  or  Y) 

0  -X", 

1  -Y* 

(float)  Detector  pitch  in  TDI  direction  in  microns 
d  >=  0.0 

(float)  the  value  of  focal  length  in  cm 
f  >=  0.0 

(integer)  Number  of  TDI  stages  in  use 
ntdi  >  0 

(float)  TDI  rate  match  =  TDI  rate  /  LOS  motion  rate 
beta  >=0.0 

(integer)  Number  of  clock  phases 
phase  >  0 


-d 

-f 

-ntdi 

-beta 

-phase 
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D.  Optics  Man  Pages 


D.l.  opdefootf  —  OTF  due  to  Gaussian  optical  blur  function 


Object  Information 

Category:  IBSM  Subcategory:  Optics 

Operator:  Defocus  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  the  OTF  due  to  optics  defocus  to  an  input  spatial  spectrum.  The 
defocus  (or  blur)  function  is  modeled  as  a  gaussian  distribution  for  which  the  full-width,  half  maximum 
(FWHM)  in  angular  space  is  specified  independently  for  the  X  and  Y  axes.  Spatial  scaling  (rad-1  units) 
and  band  attributes  must  be  available  for  the  input  image,  and  the  data  must  correspond  to  a  spatial 
spectrum,  which  can  be  produced  from  an  image  using  the  "Transform"  tool. 


Command  Line  Arguments 

Usage  for  opdefootf:  OTF  due  to  Gaussian  optical  blur  function 
%  opdefootf 

-i  (infile)  Name  of  input  IBSM  spatial  spectriim  file 

-o  (outfile)  Name  of  output  IBSM  spatial  spectrum  file 

-wx  (float)  FWHM  X  width  of  blur  function  in  mrad 
wx  >=  0.0 

-wy  (float)  FWHM  Y  width  of  blur  function  in  mrad 
wy  >=0.0 


D.2.  opdiffotf  —  OTF  due  to  aperture  diffraction 


Object  Information 

Category:  IBSM  Subcategory:  Optics 

Operator  Diffraction  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  a  diffraction  limited  OTF  to  an  input  spatial  spectrum  for  a  circular 
or  rectangular  aperture.  The  relative  obscuration  applies  only  to  a  centrally  obscured  circular  aperture, 
and  is  the  ratio  of  obscuration  to  aperture  diameters.  Spatial  scaling  (rad-1  units)  and  band  attributes 
must  be  available  for  the  input  image,  and  the  data  must  correspond  to  a  spatial  spectmm,  which  can  be 
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produced  from  an  image  using  the  "Transform"  tool. 


Command  Line  Arguments 

Usage  for  opdiffotf:  OTF  due  to  aperture  diffraction 
%  opdiffotf 

-i  (infile)  Name  of  input  IBSM  spatial  spectrum  file 

-o  (outfile)  Name  of  output  IBSM  spatial  spectrum  file 

-shape  (cycle)  Aperture  shape  description 
0  "Circular", 

1  "Rectangular" 

-dx  (float)  Aperture  diameter  (circular)  or  X  width  (rectcingular )  in  cm 

dx  >=  0.0 

-dy  (float)  Aperture  Y  width  (rectangular  only)  in  cm 

dy  >=  0.0 

-obs  (float)  Ratio  of  obscuration  to  aperture  diameter  (circular  only) 
unbounded 


D.3.  opdrifotf —  OTF  due  to  linear  line-of-sight  drift 


Object  Information 

Category:  IBSM  Subcategory:  Optics 

Operator:  Drift  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  OTF  due  to  a  linear  angular  drift  of  the  sensor  line-of-sight  dur¬ 
ing  the  integration  time  to  an  input  spatial  spectrum.  The  magnitude  of  the  linear  drift  is  specified  inde¬ 
pendently  for  the  X  and  y  directions.  Only  the  effect  of  drift  on  the  OTF  is  modeled;  possible  image  dis¬ 
tortion  effects  are  not.  Spatial  scaling  (rad-1  units)  and  band  attributes  must  be  available  for  the  input 
image,  and  the  data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced  from  an  image  using 
the  "Transform"  tool. 


Command  Line  Arguments 

Usage  for  opdrifotf:  OTF  due  to  linear  line-of-sight  drift 
%  opdrifotf 


-i 

(infile) 

Name  of  input  IBSM  spatial  spectrum  file 

-o 

(outfile) 

Name  of  output  IBSM  spatial  spectrum  file 

-ax 

(float) 

Angular  X  drift  in  mrad 

ax  >=  0. 

0 

-ay 

(float) 

Angular  Y  drift  in  mrad 

ay  >=  0.0 
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D.4.  opjittotf  —  OTF  due  to  high  frequency  line-of-sight  jitter 


Object  Information 

Category:  IBSM  Subcategory:  Optics 

Operator:  Jitter  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  OTF  due  to  high  frequency  (with  respect  to  the  sensor  integra¬ 
tion  time)  jitter  of  an  optical  sensor  line-of-sight.  By  virtue  of  this  high  bandwidth  assumption,  the  jit¬ 
ter  is  specified  by  its  x  and  y  axis  root-mean  square  (RMS)  components  and  is  modeled  with  a  zero- 
mean,  gaussian  assumption.  Spatial  scaling  (rad-1  units)  and  band  attributes  must  be  available  for  the 
input  image,  and  the  data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced  from  an  image 
using  the  "Transform"  tool. 


Command  Line  Arguments 

Usage  for  opjittotf:  OTF  due  to  high  frequency  line-of-sight  jitter 
%  opjittotf 

(infile)  Name  of  input  IBSM  spatial  spectrum  file 
Name  of  output  IBSM  spatial  spectrum  file 
RMS  X  jitter  in  mrad 


“i 

(infile) 

Name 

-o 

(outfile) 

Name 

-sigx 

(float) 

RMS  : 

sigx  >= 

0.0 

-sigy 

(float) 

RMS  ' 

sigy  >= 

0.0 

RMS  Y  jitter  in  mrad 


D.5.  opradiom  —  Conversion  of  pupil  plane  radiance  to  focal  plane  irradiance 


Object  Information 

Category:  IBSM  Subcategory:  Optics 

Operator  Radiometry  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  performs  a  radiometric  conversion  of  an  eflfective  pupil  plane  radiance  image  to  a  focal 
plane  irradiance  image.  The  model  accommodates  sensor  throughput,  thermal  radiation  of  the  optical 
chain,  and  stray  radiance.  The  input  data  must  be  properly  scaled  in  radiance  (W/m2sr)  units. 


Command  Line  Arguments 

Usage  for  opradiom:  Conversion  of  pupil  plane  radiance  to  focal  plane  irradiance 
%  Opradiom 
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-i 

(infile) 

Name  of  input  IBSM  radiance  image  file 

-o 

(outfile) 

Name  of  output  IBSM  irradiance  image  file 

“fnum 

(float) 

Effective  sensor  F -number 

fnum  >  0 

.0 

-trans 

(float) 

Mean  optical  transmission  of  sensor  over  band 

unbounded 

-em 

(float) 

Effective  emissivity  of  sensor  optics  over  bcuid 

unbounded 

-temp 

(float) 

Effective  temperature  of  sensor  optics  in  K 

temp  >= 

0.0 

-stray 

(float) 

Effective  sensor  stray  radiance  in  uW/cm2sr 

stray  >= 

0.0 

D.6.  opuserotf  —  OTF  defined  by  1-D  user  data 


Object  Infonnation 

Category:  IBSM  Subcategory:  Optics 

Operator:  User  OTF  (1-D)  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  OTF  based  on  user  data  to  an  input  spatial  spectrum.  The  data  is 
contained  in  an  ASCII  file  in  formats  accommodating  real  or  complex  data  specified  with  radial  form 
(radial),  with  separability  in  the  x  and  y  axes  (x  or  y),  or  with  separability  and  similarity  in  the  x  and  y 
axes  (x  and  y).  The  user  data  is  interpolated  using  linear  or  spline  techniques  to  the  input  spatial  spec¬ 
trum  grid  prior  to  forming  the  2-D  OTF.  Spatial  scaling  (iad-1  units)  and  band  attributes  must  be  avail¬ 
able  for  the  input  image,  and  the  data  must  correspond  to  a  spatial  spectrum,  which  can  be  produced 
from  an  image  using  the  "Transform"  tool.  The  input  file  contains  a  series  of  lines,  one  for  each  OTF 
point,  in  the  following  format: 

Radial,real:  Frequency  OTF 

Radial, complex:  Frequency  OTFreal  OTFimag 

X  or  Yo’eal:  Frequency  OTF 

X  or  Y,complex:  Frequency  OTFreal  OTFimag 

X  and  Y,real:  FrequencyX  OTFX  FrequencyY  OTFY 

X  and  Y,complex:  FrequencyX  OTFXreal  OTFXimag  FrequencyY  OTFY real  OTFTrtmag 


Command  Line  Arguments 

Usage  for  opuserotf:  OTF  defined  by  1-D  user  data 
%  opuserotf 

-i  (infile)  Name  of  input  IBSM  spatial  spectrum  file 
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-o  (outfile)  Name  of  output  IBSM  spatial  spectriim  file 

-fname  (infile)  Name  of  ASCII  file  containing  user  OTF  data 

-typel  (cycle)  User  OTF  data  format 

0  'Radial", 

1  'X  or  Y', 

2  'X  and  Y" 

-type2  (cycle)  User  OTF  data  type 

0  "Real", 

1  'Complex* 

-type3  (cycle)  Interpolation  type 

0  "Linear", 

1  "Spline" 


D.7.  opusrot£2  —  OTF  defined  by  2-D  user  data  in  OTF  or  PSF  form 


Object  Information 

Category:  IBSM  Subcategory:  Optics 

Operator:  User  OTF  {2-D)  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  OTF  based  on  bilinear  interpolation  of  a  user  2-D  OTF  or  dis¬ 
crete  Fourier  transformation  of  a  user  2-D  PSF.  The  output  data  format  (e.g.  real  vs.  complex)  is  defined 
by  the  input  spatial  spectrum  file,  which  must  contain  spatial  scaling  (rad-1  units)  attributes.  Normaliza¬ 
tion  produces  a  unity  OTF  at  zero  frequency. 


Command  Line  Arguments 


Usage  for  opusrotf2 :  OTF  defined  by  2-D  user  data  in  OTF  or  PSF  form 
%  opusrot£2 


-i 

(infile) 

Name  of  input  : 

-o 

(outfile) 

Name  of  output 

-otf 

(infile) 

Name  of  ASCII 

-type 

(cycle) 

0  -PSF-, 

1  -OTF" 

File  type 

-norm 

(cycle) 

0  "No", 

1  -Yes" 

Normalization 

-bands 

(cycle)  OTF /PSF  bands 

0  "Single", 

1  "Multiple" 
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D.8.  opwaveotf  —  OTF  due  to  random  pupil  plane  wavefront  error 


Object  Information 

Category:  IBSM  Subcategory:  Optics 

Operator:  Wavefront  OTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  and  applies  an  OTF  due  to  a  random  wavefront  error  specified  in  the  sensor 
pupil  plane  to  an  input  spatial  spectrum.  The  statistics  of  the  wavefront  error  are  modeled  as  zero  mean, 
gaussian  with  independently  specified  correlation  lengths  for  the  x  and  y  axes  in  the  pupil  plane.  Spa¬ 
tial  scaling  (rad-1  units)  and  band  attributes  must  be  available  for  the  input  image,  and  the  data  must 
correspond  to  a  spatial  spectrum,  which  can  be  produced  from  an  image  using  the  "Transform"  tool. 

Command  Line  Arguments 

Usage  for  opwaveotf:  OTF  due  to  random  pupil  plane  wavefront  error 
%  opwaveotf 

-i  (infile)  Name  of  input  IBSM  spatial  spectrum  file 

-o  (outfile)  Name  of  output  IBSM  spatial  spectrum  file 

-sigma  (float)  RMS  wavefront  error  in  waves 

sigma  >=  0.0 

-wref  (float)  Reference  wavelength  in  microns 
wref  >=  0.0 

-lx  (float)  Wavefront  error  X  correlation  length  in  cm 

lx  >=  0.0 

-ly  (float)  Wavefront  error  Y  correlation  length  in  cm 

ly  >=  0.0 


Object  Information 

Category:  IBSM  Subcategory:  Performance 

Operator  Create  MTF  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  a  ID  MTF  based  on  a  real- valued  2D  spatial  spectrum  (MTF)  file.  The 
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generated  MTF  can  lie  along  either  the  X  or  Y  axes  of  the  input  data,  and  is  placed  as  a  single  row  in  an 
IBSM  output  file.  The  frequency  scale  (min/max)  can  be  either  the  same  as  the  input  (automatic)  or 
overriden  (fixed).  In  the  latter  case ,  the  data  is  linearly  interpolated  to  the  specified  grid.  This  generated 
IBSM  file  is  of  a  form  suitable  for  input  to  the  "Sensor  Performance"  function. 


Command  Line  Arguments 


Usage  for  pecremtf:  Generate  1~D  MTF  from  2-D  MTF 
%  pecremtf 

Name  of  input  IBSM  spatial  spectmam  file 
Name  of  output  IBSM  spatial  spectziim  fil€ 
Designator  of  MTF  axis 


Designator  for  spatial  frequency  scaling  method 


Minimum  spatial  frequency  (fixed  scale) 
0.0,  default  =  0) 

Maximum  spatial  frequency  (fixed  scale) 
0.0,  default  =  100000) 

Number  of  spatial  frequency  points  (fixt 
0,  default  =  100) 


-i 

(infile) 

-o 

(outfile) 

-axis 

(cycle) 

0  “X" , 

1  -Y- 

-scale 

(cycle) 

0  "Auto", 
1  "Fixed" 

[-fmin] 

(float) 
(fmin  >= 

[-fmax] 

(float) 
(fmax  >= 

[-nfreq] 

( integer ) 
(nfreq  > 

E.2.  pecrespec  —  Generation  of  atmospheric  and  sensor  spectral  distributions 


Object  Information 

Category:  IBSM  Subcategory:  Performance 

Operator:  Create  Spectra  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  an  eight-band  IBSM  file  containing  atmospheric  and  sensor  spectral  distribu¬ 
tions.  The  bands  contain,  respectively,  (0)  atmospheric  path  transmission,  (1)  atmospheric  path  radi¬ 
ance,  (2)  irradiance  at  the  target  fro  a  direct  source  (sun  or  moon),  (3)  diffuse  downwelling  radiance  at 
the  target,  (4)  target  reflectance,  (5)  background  reflectance,  (6)  sensor  optical  chain  transmission,  and 
(7)  detector  quantum  efficiency.  All  are  in  single  row  format  of  "nwavl"  columns.  When  the  override 
parameters  are  not  selected,  the  first  four  distributions  (0-3)  are  computed  via  MODTRAN  in  reference 
to  the  ASCn  path  data  file,  the  following  two  (4-5)  from  the  ASCH  reflectance  file  data,  and  the  final 
two  (6-7)  from  the  ASCII  sensor  file  data.  In  all  cases,  the  spectral  distributions  are  computed  through 
a  linear  interpolation  of  the  data,  and  can  be  overridden  with  the  spectrally  flat  override  parameters. 
The  minimum  and  maximum  x  parameters  in  the  EBSM  header  contain  the  wavelength  minimum  and 
maximum  in  meters.  The  minimum  y,  maximum  y,  center  wavelengths,  and  spectral  bandwidths  con¬ 
tained  in  the  header  should  be  ignored.  This  generated  IBSM  file  is  of  a  form  suitable  for  input  to  the 
"Sensor  Performance"  function. 
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Command  Line  Arguments 


for  pccrespec I  Generation  of  atinospheric  and  sensor  spectral  distributions 
%  pecrespec 


-o 

-wmin 


-wmax 


-nwavl 


(out file)  Name  of  output  IBSM  spectral  distribution  file 
(float)  Minimum  wavelength  of  distributions  in  microns 
unbounded 

(float)  Maximxim  wavelength  of  distributions  in  microns 
unbounded 

(integer)  Number  of  spectral  points  in  the  distributions 


nwavl  >  0 


[-i] 

( “wor  kdir  ] 

[-fname_atm] 

(“transl 

[-pathrad] 

[-direct] 

[-zenith] 

[-diffuse] 

[-fncune_ref  1] 

[-nrefl] 

[-rtgt] 

[-rbkg] 

[  -  f  name_sens  ] 

[-nsens] 

[-topt] 

[-eta] 


(infile)  trigger  input  to  sync  execution 
(default  =  {none}) 

(string)  Optional  working  directory  for  MODTRAN 
(default  =  {none}) 

(infile)  Name  of  the  input  MODTRAN  path  data  file 
(default  =  {none}) 

(float)  Atmospheric  path  transmission  (override) 

(unbounded,  default  =  1) 

(float)  Atmospheric  path  radiance  in  uW/cm2umsr  (override) 

(pathrad  >=  0.0,  default  =0) 

(float)  Direct  Source  irradiance  at  the  target  in  uW/cin2um  (override) 
(direct  >=  0.0,  default  =  0) 

(float)  Source  zenith  angle  in  degrees  (used  with  direct  irradiance) 
(unbounded,  default  =  0) 

(float)  Diffuse  downwelling  radiance  in  uW/cm2umsr  (override) 

(diffuse  >=  0.0,  default  =  0) 

(infile)  Name  of  ASCII  reflectamce  data  file 
(default  =  {none}) 

(integer)  Number  of  lines  in  the  reflectance  data  file 
(2  <=  nrefl  <=  2.14748e+09,  default  =  2) 

(float)  Target  reflectance  (override) 

(tinbounded,  default  =  1) 

(float)  Backgroimd  reflectance  (override) 

(unbounded,  default  =  0) 

(infile)  Name  of  the  ASCII  sensor  data  file 
(default  =  {none}) 

(integer)  Nximber  of  lines  in  the  sensor  data  file 
(2  <=  nsens  <=  2.14748e+09,  default  =  2) 

(float)  Optical  chain  transmission  (override) 

(unbounded,  default  =  1) 

(float)  Detector  quantum  efficiency  (override) 

(unbounded,  default  =  1) 


E.3.  peexpctl  —  Exposure  Control 


Object  Information 

Performance 

Category:  IBSM 

Subcategory: 

Operator  Exposure  Control 

In  Cantata:  Yes 

Object  Type: 

kroutine 
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Description 

This  function  performs  exposure  control  by  determining  the  pixel  dwell  time,  TDI  stages,  and  aggre¬ 
gation  levels  which  place  the  maximum  of  the  target  or  background  signal  levels  between  specified 
mnimum  and  maximum  levels.  For  the  signal  computation,  sensor  spectral  distributions  using  the  "Cre¬ 
ate  Spectra"  function  are  used.  The  output  consists  of  sensor  configuration  values  which  are  placed  in  a 
5  band  floating-point  VIFF  file  (one  value  per  band),  which  canbe  extracted  to  global  variables  using 
the  "Print  Value"  function.  The  output  file  format  is  as  follows: 


Band  0  Pixel  dwell  time  in  msec 

Band  1  TDI  stages 

Band  2  X  aggregation 

Band  3  Y  aggregation 

Band  4  Signal  level  (max)  in  electrons 


Command  Line  Arguments 


Usage  for  peexpctl :  Exposure  Control 
%  peexpctl 

_i_spect  (infile)  Name  of  input  IBSM  spectral  distribution  file 
_o  (outfile)  Name  of  output  VIFF  result  file 

-dwell  (float)  Desired  pixel  dwell  time  in  msec 

dwell  >=  0.0 

-tdi  (string)  Allowable  TDI  stages  (input  as  string  with  blanlc  separators) 

-agx  (string)  Allowable  X  aggregations  (input  as  string  with  blank  separators) 

-agy  (string)  Allowable  Y  aggregation  (input  as  string  with  blank  separators) 

-range  (float)  Target  slant  range  in  km 

range  >  0.0 

-ttgt  (float)  Target  tenperature  in  Kelvin 

ttgt  >=  0.0 

-tbkg  (float)  Background  tenperature  in  Kelvin 

tbkg  >=  0.0 

-f  (float)  Sensor  focal  length  in  cm 

f  >=  0.0 

•  -fnxim  (float)  Effective  sensor  F-number 

fnum  >=  0.0 

-eopt  (float)  Effective  emissivity  of  sensor  optics  over  band 

unbounded 

-topt  (float)  Effective  tenperature  of  sensor  optics  in  Kelvin 

topt  >~  0.0 

-Istray  (float)  Effective  sensor  stray  radiance  in  uW/cm2sr 

Istray  >=  0.0 

-voc  (float)  Detector  X  width  in  microns 

wx  >=  0.0 

(float)  Detector  Y  width  in  microns 
wy  >=  0.0 

-jd  (float)  Detector  dark  current  density  in  microaB5>s/cm2 

jd  >=  0.0 

-nmin  (float)  Minimum  signal  level  for  exposure  control  in  electrons 

nmin  >=0.0 

-nmax  (float)  Maucimum  signal  level  for  exposure  control  in  electrons 

nmax  >0.0 


-agx 

-agy 

-range 

-ttgt 

-tbkg 


-fnum 


-eopt 


-Istray 
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E.4.  peimgqual  —  Computes  edge  response  functions  and  NIIRSfrom  test  image 


Object  Information 

Category:  IBSM  Subcategory:  Performance 

Operator:  Image  Quality  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  computes  image  quality  metrics  including  edge  response  functions,  GSD,  SNR,  and 
NIIRS  from  an  image  containing  either  test  patterns  or  edge  features.  The  user  specifies  the  locations  of 
horizontal  and  vertical  edges  as  well  as  uniform  image  regions  for  noise  estimation.  Outputs  include 
ASCII  files  containing  the  X  and  Y  edge  response  functions  (Two  column  format)  and/or  a  metric  sum¬ 
mary  file. 


Command  Line  Arguments 


for  peiiHQQUdl :  Con5)ut.es  ed^e  response  functions  and  NIIRS  from  test  linage 


%  peimgqual 
-i 

~xve 


-yve 


-hve 


“Xhe 


-yhe 


-hhe 


-xtr 

-ytr 


-xbr 


-ybr 


-vxs 


-hrs 

-range 


(infile)  Input  IBSM  image  file  name 
(float)  X  position  of  vertical  edge  in  mrad 
xanbounded 

(float)  Y  position  of  vertical  edge  in  mrad 
unbounded 

(float)  Vertical  edge  size  in  mrad 
unbounded 

(float)  X  position  of  horizontal  edge  in  mrad 
\anbounded 

(float)  Y  position  of  horizontal  edge  in  mrad 
unbounded 

(float)  Horizontal  edge  size  in  mrad 
unbounded 

(float)  X  position  of  target  region  in  mrad 
unbounded 

(float)  Y  position  of  target  region  in  mrad 
\inbounded 

(float)  X  position  of  background  region  in  mrad 
unbounded 

(float)  Y  position  of  background  region  in  mrad 
unbounded 

(float)  Vertical  region  size  in  mrad 
unbounded 

(float)  Horizontal  region  size  in  mrad 
unbounded 

(float)  Nominal  slant  rauige  in  km 
range  >  0.0 


[-ol]  (out file)  Output  ASCII  X  response  file  name 
(default  =  {none}) 

[-02]  (outfile)  Output  ASCII  Y  response  file  name 
(default  =  (none) ) 

[-o3]  (outfile)  Output  ASCII  results  file  name 
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(default  =  {none}) 


E.5.  pesensper  —  Compute  sensor  performance  metrics 


Object  Information 

Category:  IBSM  Subcategory:  Performance 

Operator:  Sensor  Performance  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  computes  metrics  characterizing  sensor  performance  (including  SNR,  GRD,  and  NEKS) 
based  on  a  set  of  sensor  parameters;  target,  background,  atmosphere,  and  sensor  spectral  distributions 
using  the  "Create  Spectra"  function;  and  x  and  y  axis  MTF  functions  using  the  "Create  MTF"  function. 
The  metric  is  written  to  an  ASCII  file  along  with  an  arbitrary  independent  variable.  Depending  on  the 
’update’  parameter,  the  file  is  initialized  before  writing  or  the  data  is  appended  to  an  existing  file.  In  the 
long  form  output  data  format,  the  full  input  parameter  set  is  written  to  the  file  upon  initialization 
(update  =  0),  and  the  metric  data  is  augmented  with  additional  intermediate  computation  results.  The 
short  form  output  data  format  results  in  a  two  column  ASCII  file  containing  the  independent  variable 
and  computed  metric  for  each  update  execution,  which  is  useful  for  generating  parametric  plot  data. 
The  SNR  metric  refers  to  the  target  to  background  difference  signal  relative  to  the  detected  noise  level. 
In  this  case,  the  axis  is  ignored.  The  GRD  metric  is  provided  in  meters  and  is  given  in  a  plane  orthogo¬ 
nal  to  the  sensor  line-of-sight  axis.  The  image  quality  metric  (NIIRS)  is  also  computed  relative  to  an 
orthogonal  plane  ground  sample  distance  (GSD). 


Command  Line  Arguments 


Usage  for  pesensper:  Compute  sensor  performance  metrics 


%  pesensper 
-inspect 
-o 

-metric 

0 

1 

2 


(infile)  Name  of  input  IBSM  spectral  distribution  file 
(out file)  Name  of  output  ASCII  result  file 
(cycle)  Sensor  perfonn2mce  metric 
"Signal  to  Noise  Ratio", 

"Ground  Resolved  Disteuice", 

"Image  Quality" 


-wx 

-xvar 

-wy 

-update 


-px 


-format 


0 

1 


(float)  Detector  X  width  in  microns 
wx  >=  0.0 

(float)  Independent  variaOble  for  plots 
unbounded 

(float)  Detector  Y  width  in  microns 
wy  >=  0.0 

(integer)  Output  file  update  switch  (0  =  new  file,  1  =  update) 
0  <=  update  <=  1 

(float)  Detector  X  pitch  in  microns 
px  >=  0.0 

(cycle)  Output  data  format 

"Short  Form", 

"Long  Form* 
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-py 

-td 

-rauige 

-jd 

-ttgt 

-ntdi 

-tbkg 

-agx 

-agy 

-nmin 

-f 

-nmax 

-fnum 

-eopt 

-bits 

-topt 

-sigma 

-Is tray 

-gain 


(float)  Detector  Y  pitch  in  microns 
py  >=  0.0 

(float)  Detector  dwell  (integration)  time  in  msec 
td  >=  0.0 

(float)  Target  slant  range  in  )cm 
range  >  0.0 

(float)  Detector  dark  current  density  in  inicroairps/cm2 
jd  >=  0.0 

(float)  Target  tenperature  in  Kelvin 
ttgt  >=0.0 

(integer)  Nvimber  of  detector  TDI  stages 
ntdi  >  0 

(float)  Backgroiand  temperature  in  Kelvin 
tbkg  >=  0.0 

(integer)  Aggregation  in  X  direction 
agx  >  0 

(integer)  Aggregation  in  Y  direction 
agy  >  0 

(float)  Qucuitizer  minimum  level  in  electrons 
nmin  >=  0.0 

(float)  Sensor  focal  length  in  cm 
f  >=  0.0 

(float)  Quantizer  maximum  level  in  electrons 
nnicuc  >  0.0 

(float)  Effective  sensor  F-number 
fnum  >=  0.0 

(float)  Effective  emissivity  of  sensor  optics  over  band 
unbounded 

(integer)  Digitization  levels  as  power  of  2 
bits  >  0 

(float)  Effective  tenperature  of  sensor  optics  in  Kelvin 
topt  >=  0.0 

(float)  Additive  noise  level  RMS  in  electrons 
sigma  >=  0,0 

(float)  Effective  sensor  stray  radiance  in  uW/cm2sr 
Istray  >=  0.0 

(float)  Processing  noise  gain 
gain  >  0.0 


I-i_mtfx]  (infile)  Name  of  input  IBSM  MTF  file  (X  eucis) 

(default  =  {none)) 

[-i_mtfy]  (infile)  Name  of  input  IBSM  MTF  file  (Y  axis) 

(default  =  (none)) 

[-age]  (float)  Automatic  gain  correction  level  in  electrons  (optional) 

(age  >  0.0,  default  =  le+07) 


F.  Processing  Man  Pages 
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F.l.  printerlace  —  Interlaces  IBSM  images 


Object  Information 

Category:  IBSM  Subcategory:  Processing 

Operator:  Image  Interlace  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  interlaces  up  to  four  images  in  the  user  specified  direction:  X  (width),  Y  (height),  or  band 
(element).  The  interlaced  images  must  all  be  the  same  size. 


Command  Line  Arguments 

Usage  for  printerlace:  Interlaces  IBSM  images 
%  printerlace 


-il 

(infile) 

First  IBSM  Image  to  Interlace 

(Required) 

-i2 

( infile) 

Second  IBSM  Image  to  Interlace 

(Required) 

-o 

(outfile) 

Output  IBSM  Interlaced  File 

-mode 

0 

(cycle) 

-  X 

Interlace  Mode  Selection 

1 

-  Y-, 

2 

•’Bands" 

[-13] 

(infile) 

Third  IBSM  Image  to  Interlace 

(Optional) 

(default 

=  {none} ) 

[-14] 

(infile) 

Fourth  IBSM  Image  to  Interlace 

(Optional) 

(default 

-  {none}) 

F.2.  praggreg  —  Sums  image  elements  in  width,  height,  and  bands 


Object  Information 

Category:  IBSM  Subcategory:  Processing 

Operator  Aggregate  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  sums  image  elements  in  width,  height,  and  bands  to  a  user  specified  aggregation  level.  By 
virtue  of  the  aggregation,  the  output  image  is  smaller  in  each  dimension  by  the  aggregation  factor.  The 
residual  image,  when  the  size  is  not  an  integer  multiple  of  the  aggregation  factor,  is  truncated.  If  the 
projection  option  is  selected  for  a  particular  dimension,  the  aggregation  factor  is  overriden  by  the  by  the 
dimension  size.  Note  that  overflows  may  occur  in  the  summation  process. 


128 


Man  Pages 


IBSM  Manual  -  Chapter 


Command  Line  Arguments 


Usage  for  praggreg;  Sums  image  elements  in  width,  height,  and  beinds 


%  praggreg 
-i 
-o 

-xagg 

-yagg 

-bagg 

“projx 

-projy 

-projb 


(infile)  IBSM  file  used  as  input  to  the  aggregation  function 
(outfile)  Output  IBSM  residual  image  from  the  aggregation  function 
(integer)  User  specified  aggregation  factor  in  X 
1  <=  xagg  <=  2.14748e+09 

(integer)  User  specified  aggregation  factor  in  Y 
1  <=  yagg  <=  2.14748e+09 

(integer)  User  specified  aggregation  factor  in  bcuids 
1  <=  bagg  <=  2.14748e+09 

(boolean)  Project  input  image  in  X  (Overrides  X  Aggregation  Factor  when  true) 

(boolean)  Project  input  image  in  Y  (Overrides  Y  Aggregation  Factor  when  true) 

(boolean)  Project  input  image  in  bands  (Overrides  Band  Aggregation  Factor  when 


G.  Utilities  Man  Pages 


G.l.  ibsminfo  —  Extract  IBSM  header  information 


Object  Infonnation 

Category:  IBSM  Subcategory:  Utilities 

Operator:  List  IBSM  File  Data  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  generates  an  ASCII  file  containing  infonnation  extracted  from  the  IBSM  header  and 
attributes.  This  includes  data  type,  image  size,  number  of  bands,  spatial  scaling,  and  the  center 
wavelength  and  spectral  bandwidth  for  each  band. 


Command  Line  Arguments 

Usage  for  ibsminfo:  Extract  IBSM  header  infonnation 
%  ibsminfo 

-i  (infile)  Name  of  input  IBSM  file 

[-o]  (outfile)  Name  of  output  ASCII  file  (optional) 
(default  =  {none}) 
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G.2.  utcreplot  —  Formats  IBSM  Data  for  Plotting 


Object  Information 

Category;  IBSM  Subcategory:  Utilities 

Operator:  Create  Plot  Data  Object  Type:  kroutine 

In  Cantata;  ^es 

Description 

This  function  extracts  a  row  or  column  of  an  IBSM  file  and  outputs  it  along  with  spatial  scaling  data 
into  two  column  ASCII  format  for  plotting. 


Command  Line  Arguments 


Usage  for  utcreplot:  Formats  IBSM  Data  for  Plotting 
%  utcreplot 


(infile) 

input  IBSM  image  or  spectrum 

-o 

(outfile) 

resulting  output  ASCII  plot  file 

-band 

(integer) 

the  value  of  integer 

0  <=  band  <=  1023 

-direct 

(cycle) 

image  direction  cycle  value 

0  "Row" , 

1  "Column" 

-rcnum 

( integer) 

the  value  of  row  or  column  number 

0  <=  rcnum  <=  32768 


G3.  utcreradi  —  Convert  Reflectance/Temperature  Image  to  Radiance  Image 


Object  Information 

Category:  IBSM  Subcategory:  Utilities 

Operator:  Create  Radiance  Image  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  converts  a  ground  reflectance  and,  optionally,  temperature  distribution  into  a  radiance 
image,  including  both  reflected  and  radiated  components.  The  reflected  component  uses  only  the 
reflectance  image  input  (assumed  Lambertian)  with  direct  and  diffuse  illumination  computed  via  MOD- 
TRAN  in  reference  to  the  ASCII  path  data  file.  This  file  can  be  generated  using  the  "MODTRAN 
Input"  tool.  It  is  important  that  the  wavelength  limits  specified  in  the  MODTRAN  path  data  file  cover 
each  spectral  band  specified  in  the  input  image  header.  The  direct,  diffuse,  or  both  illumination  values 
from  MODTRAN  can  be  overridden.  The  radiated  component  is  computed  using  both  the  reflectance 
and  temperature  (Kelvin  units)  images,  and  is  ignored  when  the  latter  does  not  exist.  The  output  image 
has  radiance  units  (W/m3sr). 
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Command  Line  Arguments 


Usage  for  utcreradi:  Convert  Ref lec tame e /Temperature  Image  to  Radiance  Image 
%  utcreradi 

-i_refl  (infile)  Name  of  input  IBSM  reflectance  image  file 

-o_radi  (outfile)  Name  of  output  IBSM  radiance  image  file 


[-i^temp] 

[-workdir] 

[-ijmod] 

[-direct] 

[-zenith] 

[-diffuse] 


(infile)  Name  of  input  IBSM  tenperature  image  file  (optional) 
(default  =  {none}) 

(string)  Optional  working  directory  for  MODTRAN 
(default  =  (none)) 

(infile)  Name  of  input  ASCII  MODTRAN  path  data  file 
(default  =  {none}) 

(float)  Direct  source  irradiance  at  target  in  uW/cm2um  (override) 
(direct  >=  0.0,  default  =  0) 

(float)  Source  zenith  angle  in  degrees  (used  with  the  direct  field) 
(unbounded,  default  =  0) 

(float)  Diffuse  downwelling  radiance  in  uW/cm2umsr  (override) 
(diffuse  >=  0.0,  default  =  0) 


G.4.  utgrdangl  —  Transform  of  image  between  ground  and  line-of-sight  angle  domains 


Object  Information 

Category;  IBSM  Subcategory:  Utilities 

Operator:  Ground/Angle  Transform  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  performs  a  coordinate  transformation  of  a  radiance  image  between  the  horizontal  ground 
plane  (meter  units)  and  line-of-sight  angular  (radian  units)  coordinates.  Only  the  image 
scaling  attribute  values  are  altered  based  on  the  viewing  geometry.  The  defined  azimuth  axis  is  the 
non-  squinted  axis  (normal  to  the  line-of-sight  and  vertical  axes).  The  scaling  of  the  squinted  axis  is 
based  on  a  linear  projection,  which  occurs  under  the  assumption  that  the  ground  field-of-view  in  the 
squinted  direction  is  much  smaller  than  the  slant  range.  No  image  distortion  or  perspective  effects  are 
modeled. 


Command  Line  Arguments 


Usage  for  utgrdangl;  Transform  of  image  between  ground  and  line-of-sight  angle  domains 
%  utgrdangl 

-i  (infile)  Name  of  input  IBSM  radiance  image  file 

-o  (outfile)  Name  of  output  IBSM  radiance  image  file 

-range  (float)  Target  slant  range  in  km 


-zenith 

-azimuth 

0 

1 


remge  >=0.0 

(float)  Target  zenith  angle  in  degrees  (downlooking  =  180) 
unbounded 

(cycle)  Definition  of  azimuth  (non- squinted)  axis 

"X-, 

•Y- 
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-direction  (cycle)  Transform  direction 
0  “Ground  to  Angle", 

1  "Angle  to  Ground" 


G.5.  utloadvars  —  Load  Cantata  Variables  from  a  File 


Object  Information 

Category:  IBSM  Subcategory:  Utilities 

Operator:  Load  Variables  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  loads  a  set  of  global  variable  definitions  from  a  user  specified  ASCII  file.  The  file  format 
parallels  the  global  variables  input  form,  and  can  include  defining  variables  as  constants,  expressions, 
and  comments  (starting  with  #). 


Command  Line  Arguments 

Usage  for  utloadvars:  Load  Cantata  VarieUoles  frcxn  a  File 
%  utloadvars 

-i  (infile)  Loads  the  Cantata  variables  from  the  input  file 


G.6.  utxfomi  —  Transforms  IBSM  Files  to/from  Image/Spectrum 


Object  Information 

Category:  IBSM  Subcategory:  Utilities 

Operator.  IBSM  Transform  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  performs  conversions  between  the  image  and  spatial  spectram  domains  via  2-D  Fast 
Fourier  Transform  (FFT)  techniques.  In  addition  to  the  FFT,  the  function  performs  scale  factor  conver¬ 
sion  of  the  header  quadrant  swapping  in  both  domains,  and  zero  padding  to  the  transform  size. 
The  transform  size  must  be  a  power  of  two  in  both  directions,  and  must  be  greater  than  or  equal  to  the 
corresponding  input  image  or  spatial  spectrum  size. 


Command  Line  Arguments 


Usage  for  utxfonn:  Transforms  IBSM  Files  to/from  Image /Spectrum 
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%  utxfonn 
-i 
-o 

-nnx 

-nny 

-otype 

1 

2 

3 

4 

5 


(infile)  Name  of  input  IBSM  image  or  spatial  spectrum  file 
(out file)  Name  of  output  IBSM  spatial  spectrum  or  image  file 
(integer)  Transform  size  (saitples)  in  X  direction 
1  <=  nnx  <=  16384 

(integer)  Transform  size  (sait^les)  in  Y  direction 
1  <=  nny  <=  16384 
(cycle)  Output  data  type 
■Coitplex" , 

•Real", 

•Imaginary" , 

•Magnitude" , 

•Phase" 


-xform_dir  (cycle)  Trcuisform  direction 

1  "Spectrvim  to  image", 

2  "Image  to  spectrxim" 


G.7.  utxgraph  —  Generate  a  command  line  to  drive  Xgraph 


Object  Infonnation 

Category:  IBSM  Subcategory:  Utilities 

Operator  Xgraph  Plot  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  provides  a  user  interface  to  the  Xgraph  plotting  application  and  allows  up  to  four  plots  per 
graph  with  preformatted  labels  and  axis  ranges.  The  ^ta  is  input  as  separate  ASCII  files  in  two  column 
format 


Command  Line  Arguments 

Usage  for  utxgraph:  Generate  a  command  line  to  drive  Xgraph 
%  utxgraph 

[-il]  (infile)  First  Xgraph  input  file 

(default  =  {none}) 

[-Ibll]  (string)  Input  file  #1  label 
(default  =  (none)) 

[-12]  (infile)  Second  Xgraph  input  file 

(default  =  (none)) 

[-lbl2]  (string)  Input  file  #2  label 
(default  =  (none)) 

[-13]  (infile)  Third  Xgraph  input  file 

(default  =  {none}) 

[-lbl3]  (string)  Input  file  #3  label 
(default  =  {none}) 

[-14]  (infile)  Fourth  Xgraph  input  file 

(default  =  {none}) 

[-lbl4]  (string)  Input  file  #4  label 
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(default  =  {none}) 

[-15]  (infile)  Fifth  Xgraph  input  file 

(default  =  (none)) 

[-lbl5]  (string)  Input  file  #5  label 
(default  =  (none)) 

[-16]  (infile)  Sixth  Xgraph  input  file 

(default  =  {none}) 

[-lbl6]  (string)  Input  file  #6  leUoel 
(default  =  {none} ) 

[-17]  ( infile)  Seventh  Xgraph  input  file 

(default  =  {none}) 

[-lbl7]  (string)  Input  file  #7  label 
(default  =  {none}) 

[-18]  (infile)  Eighth  Xgraph  input  file 

(default  =  {none}) 

[-lbl8]  (string)  Input  file  #8  label 
(default  =  {none}) 

[-logx]  (flag)  Flag  to  select  logarithmic  X  axis 

[-logy]  (flag)  Flag  to  select  logarithmic  Y  axis 

[-xlabel]  (string)  Text  for  X  axis  label 
(default  =  {none}) 

[-ylabel]  (string)  Text  for  Y  axis  label 
(default  =  {none}) 

[-title]  (string)  Text  for  graph  title 
(default  =  {none}) 

[Optional  Mutually  Inclusive  Group  -  Specify  all  or  none  of  the  2  options:] 
[-xmin]  (float)  Minimum  X  scale  value 
(unbounded,  default  =  0) 

[-xmax]  (float)  Maximum  X  scale  value 
(unbounded,  default  =  1) 

[Optional  Mutually  Inclusive  Group  -  Specify  all  or  none  of  the  2  options:] 
[-ymin]  (float)  Minimum  X  axis  scale  value 
(unbounded,  default  =  0) 

[-ymax]  (float)  Minimum  X  axis  scale  value 
(unbounded,  default  =  1) 


G.8.  vifI2ibsm  —  Converts  a  VIFF  File  to  IBSM  Format 


Object  Infonnation 

Category:  IBSM  Subcategory:  Utilities 

Operator  VIFF  to  IBSM  Object  Type:  kroutine 

In  Cantata:  Yes 

Description 

This  function  adds  attribute  infonnation  to  a  VIFF  file  to  create  an  IBSM  format  output.  The  spatial 
scaling  and  spectral  band  data  can  either  be  entered  from  the  screen  (up  to  four  bands),  or  input  from  an 
ASCn  header  file.  The  file  format  corresponds  to  that  generated  using  the  "List  IBSM  File  Data"  func¬ 
tion. 
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Command  Line  Arguments 

Usage  for  viff2ibsm:  Converts  a  VIFF  File  to  IBSM  Format 
%  viff2ibsm 

“i  { infile)  Name  of  input  VIFF  image 

-o  (out file)  Name  of  output  IBSM  image 

-xmin  (float)  Minimum  X  value  (normally  rad  or  rad-1  units) 

tanbounded 

-xmax  (float)  Maiximum  X  value  (normally  rad  or  rad-1  units) 
unbounded 

-ymin  (float)  Minimum  Y  value  (normally  rad  or  rad-1  units) 
unbounded 

-ymax  (float)  Maximum  Y  value  (normally  rad  or  rad-1  \inits) 
unbounded 

-wll  (float)  Center  wavelength  of  band  #1  in  microns 
unbounded 

-bwl  (float)  Spectral  bandwidth  of  band  #1  in  microns 

bwl  >  0.0 

-wl2  (float)  Center  wavelength  of  band  #2  in  microns 
unbounded 

-bw2  (float)  Spectral  bandwidth  of  beuid  #2  in  microns 
bw2  >  0.0 

-wl3  (float)  Center  wavelength  of  band  #3  in  microns 
unbounded 

-bw3  (float)  Spectral  beuidwidth  of  band  #3  in  microns 
bw3  >  0.0 

-wl4  (float)  Center  wavelength  of  band  #4  in  microns 
unbounded 

-bw4  (float)  Spectral  bsuidwidth  of  band  #4  in  microns 
bw4  >  0.0 

[-iasc]  (infile)  Name  of  input  ASCII  header  file  (optional) 
(default  =  {none}) 
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